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INTRODUCTION 


(The primary object of thermodynamics is to describe the 
behaviour of systems in equihbrium, in such a way that this 
behaviour is known quantitatively when equilibrium con- 
ditions such as temperature and pressure are altered. Using 
the first and second laws of thermodynamics, it is possible 
to define certain functions, such as the entropy and free 
energy of the system, which give a formal description of the 
behaviour of the system in a very elegant way. Thermo* 
dynamics is, however, essentially a practical science. The 
:hief progress in the last thirty years has been in devising 
methods of evaluating thermodynamic functions rather than 
in extending their formal applications. From a practical 
point of view the problem is always to express the thermo- 
dynamic functions in terms of quantities which are more 
easily observed experimentally. Such progress has been 
made that the historic method of deriving the law of mass 
action from the theory of the velocity of reactions,, and 
the concept of dynamic equilibrium, has been inverted. 
One of the tentative applications of thermodynamic func- 
tions at the present time is in the calculation of reaction 
velocity. 

The aim of the present hook is to describe briefly the formal 
applications of thermodynamic functions to the calculation 
of physico-chemical equilibria. This is followed by a descrip- 
tion of the experimental methods of evaluation that are 
available. Direct measurements on equilibrium processes 
have been very fully described elsewhere (e.g. Lewis and 
Randall, Thermodynamics, McGraw Hill Co., 1923), so that 
attention is principally paid to the fundamental departure 
of evaluating entropy and free energy from the heat content, 
which is much more easily measured. This involves a careful 
discussion of the Nernst Heat Theorem. 
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viii 

Such a discussion is not possible without at least an ele- 
mentary treatment of statistical theory. The relation between 
the thermodynamic functions and the partition function 
(Zustandsumme) is therefore discussed in an elementary way, 
the aim being to give a sufficiently complete account for 
thermodynamical requirements without assuming the bur- 
dens of a more detailed statistical theory. 

The theory of the partition function is also required for 
the recently developed methods of calculating thermo- 
dynamic functions from spectroscopic data. Here again an 
attempt has been made to discuss the principal methods 
used, so as to give an account which is fairly complete from 
the thermodynamical aspect, without going too far into 
difficulties which merely arise from complications in certain 
band spectra. 

The application of velocity considerations to the calcula- 
tion of equilibria, as in the adsorption of gases, is also briefly 
described, in order to complete the list of subjects which are 
not always treated very adequately in ordinary text-books on 
thermodynamics. For the converse reason, practically no 
reference is made to the subject of electrolytes. 

The book is based on lectures delivered in the Department 
of Thermodynamics in Oxford. 



THE USE OF THERMODYNAMIC FUNCTIONS 

Thermodynamics deals primarily with, systems in equili- 
brium. It aims at describing completely the behaviour of 
all imaginable systems in equilibrium, when various factors 
such as temperature, pressure, and concentration of the 
various components are altered. With ninety-two elements 
and the innumerable compounds which can be obtained from 
them the number of imaginable systems is, however, very 
large. Even if only those of practical importance had each 
to be investigated experimentally before anything could be 
said of them, the prospect would be rather dismal. 

Fortunately, thermodynamics provides a number of theo- 
retical propositions which enormously simplify the problem. 
One of these is the Phase Rule, which states that the number 
of undetermined variables, or degrees of freedom, F , is not 
arbitrary, but is given by the rule F = G— P+2, where G 
is the minimum number of components required to describe 
each phase of the system completely, and P is the number 
of phases. 

A second way in which the study of systems in equilibrium 
can be greatly simplified is by the use of thermodynamic 
functions, which are based on the first and second laws of 
thermodynamics, and which can be evaluated for each single 
substance separately. As is explained in what follows, the 
thermodynamic functions of the substances A, B, C permit 
the evaluation of the functions for any system composed of 
them, without further experiment, so that it is possible in 
principle to calculate the position of equilibrium in all com- 
binations of A, B, G taken two or more at a time, and thus 
to cover n\ systems (apart from limitations of valency) 
with only n determinations. Although the calculation of 
thermodynamic functions offers both theoretical and experi- 
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2 USE OF THERMODYNAMIC FUNCTIONS 

mental difficulties, their very practical use, once they are 
determined, is clear. Practically all the progress in modern 
thermodynamics depends on the development of convenient 
methods for determining the thermodynamic functions of 
individual substances. 

Thermodynamic functions based on the first law 

The heat content and the internal energy E. 

Two important properties of all thermodynamic functions 
may be illustrated from one of the simplest, the heat content 
$. If two different states of a single substance or a system 
can be characterized by different values of temperature, 
pressure, concentrations, etc., the increase in heat content 
of a system in changing from state A to state B is defined as 
the energy taken up when the change is carried out at con- 
stant pressure. The symbol i and the name Enthalpy have 
also been used for this function. 

Acquisitive convention as to signs . 

It should be noted at the outset that various conventions 
as to whether a change in heat content is to be counted 
negative or positive have been used, according to the pre- 
ferences and nationality of the author. This variability leads 
to confusion. By analogy with the universal nature of human 
acquisitiveness, which transcends national boundaries, when 
a system gains in entropy, energy, etc., the change will be 
accounted as positive. In the same way mass action con- 
stants are written with the products in the numerator. The 
definition of the change in heat content may be symbolized 
by the equation A$ = — § a , where and are the heat 

contents in the states A and B. The first important property 
of the function § is that the value of A§ depends only on the 
initial and final states of the single substance or system under 
discussion, and not on the process by which the change is 
brought about. 
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This property follows from the first law of thermodynamics, 
which states that energy can neither be created nor destroyed. 
If it were possible to have two different processes with 
different changes in heat content in going from A to B , it 
would be possible by working one of the processes in the 
opposite direction to set up a cycle, leading to the indefinite 
production or destruction of energy without permanent 
change in the working substance. This is contrary to the 
first law. The value of the function § in any state is thus 
dependent only on what that state is, and not on previous 
history. 

In the language of the differential calculus, this property, 
which is common to all thermodynamic functions, though 
for different reasons, is expressed by the statement that 
d$r>,dS 3 etc., are complete differentials, i.e. that in any small 
change in the system, which can be expressed in terms of 
changes in the variables x ± , x 2 , x 3 ,... characterizing it, the 
change d§> is not arbitrary, nor dependent on previous history, 
but is given completely by the expression 



A second important property of the heat content, and of 
other thermodynamic quantities such as the volume, the 
entropy, etc., is that they are extensive, i.e. that the numerical 
value is directly proportional to the mass of working sub- 
stance taken. This follows from the fact that the change in 
heat content in any process is the same if the substance 
as a whole is brought from state A to state B, or if it is 
first subdivided (which does not require any work) and the 
individual portions are brought from state A to state B 
separately. This property may be symbolized by the equa- 
tion § = nH> where H is the heat content per mole, and there 
are n moles in the mass of working substance chosen. Exten- 
sive properties should be contrasted with intensive properties 
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such as temperature and pressure, which must be the same 
throughout a system in equilibrium and do not depend on 
the quantity of working substance. 

Hess's law . 

From these two properties it follows that the heat content 
per mole may be treated like any other quantity in a 
chemical equation. Thus from the equations 

CaC0 3 CaO+C0 2 +AH ± , 

C+C0 2 ->2C0 +A H 2 

it follows by direct addition that the change in heat content 
in the reaction 

CaC0 3 +C -> CaO+2CO +A H z 
is simply A H 3 = AH^AH^ 

Hess’s law is of great value in thermochemistry for the 
calculation of changes in heat content which cannot be 
obtained by direct calorimetry. Since all the other functions 
to be defined possess the same two properties as the heat 
content, one means of calculating thermodynamic quantities 
not directly accessible to experiment is evidently to combine 
two or more equations for which they can be est im ated. 


Measurements of heat content ; choice of a reference state. 

The change in heat content in a large number of processes 
can be measured quite conveniently in suitable forms of 
calorimeters, as well as by other means such as spectroscopic 
observations. The heat content in any single state is to some 
extent arbitrary, since only differences in heat content can 
be measured. Any one state of a substance or system may 
therefore be selected as reference state, and the heat content 
of all other states referred to the selected state as zero. 


When a suitable reference state has been chosen the change 
in heat content in passing to any other state A may be written 

in full A TT TT TT 

LXL1 a 11 a 11 Standard* 
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Where no ambiguity can arise the symbol H a is simply used 
for the heat content in the state A , referred to the standard 
state as zero. 

The particular state chosen as zero makes no difference to 
the change in heat content in proceeding from one state A to 
another B , since the difference 

^^a,b = ^standard) ifia -^standard) 

is independent of what value is chosen as standard. The 
choice of a suitable zero state for thermodynamic functions 
is, however, important when these are to be correlated with 
statistical theory, and is referred to again below. A further 
practical consideration in the choice is that the changes can 
be accurately measured, i.e. that the standard state is experi- 
mentally accessible. 

The internal energy E. 

When a system passes from state A to state B at constant 
volume the increase in internal energy is defined as 
A E = E b — E a , 

the net energy absorbed. Just as with the heat content, the 
values of E b and E a are to some extent arbitrary, and the 
internal energy of a system in any state must be referred to 
a suitable reference state, from which the differences can be 
measured. From the first law of thermodynamics, the in- 
crease in heat content in any process must be equal to the 
increase in internal energy + the work done on the system 
by the external pressure. The relation between the two 
functions can thus be written 

AH = A E+P(V b -V a ) 
or H = E+PV . 

In general, thermodynamic functions referring to changes at 
constant pressure are more convenient for experimental 
determination, whereas functions referring to constant 
volume are more suitable for theoretical calculation. The 
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relations between the different functions, which are used for 
calculating one from another, are discussed below. 

Thermodynamic functions based on the second law 

Tor calculating the equilibria among chemical substances 
it is necessary to use functions based on the second law as 



well as the first. The simplest of these is the entropy 8. The 
increase in entropy of a system in any change is defined as 
A 8 — QjT, where Q is the heat absorbed and T is the 
absolute temperature, provided the change is reversible. 
This provision is made to preserve the essential property 
that the entropy difference between any state and a standard 
reference state shall depend only on that state, and not on 
the previous history of the working substance, i.e. that dS 
shall be a complete differential. As a result of the condition 
of reversibility the quantity TAS for a change from a state 
A to a state B is in general quite different from the calori- 
metric heat of reaction, which would be observed when the 
change is carried out under irreversible conditions. 

The way in which this property depends on the second law 
of thermodynamics may be seen from a simple Carnot’s cycle. 
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AB , CD are isothermals, BC and AD adiabatics. The 
entropy change in proceeding from A to B is 

Qab 


AS 




From B to C there is no entropy change. From C to D the 
change is q 


AS = S D —S 0 


PCD 
% * 


From D to A there is no entropy change. 

From the second law of thermodynamics, if the process is 
carried out reversibly in every part, 

Qab Qcd 

T x T 2 ’ 

i.e. Sg Sjl — S D — S G . 

Provided the heat absorbed is measured in any other 
reversible change, s 

J* = S B ~S A 

A 

must still be the same. If this were not the case, it would be 
possible to construct a cycle for converting the heat of the 
surroundings indefinitely into work, contradicting the second 
law. It follows that the entropy change can depend only on 
the initial and final states, and not on the (reversible) path. 
This is the property required. 

Like the heat content and internal energy, S is an exten- 
sive property whose value is proportional to the working 
mass of the substance, so that entropy changes can be in- 
cluded as one of the terms of a chemical equation. Unlike 
the heat content, however, whenever the change has to be 
measured experimentally care must be taken that the pro- 
cess chosen is reversible at every stage. 

All the other thermodynamic functions used are derived 
from these three, and consequently share the same properties 
of depending only on the actual state of the substance (and 
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the reference state), and of being extensive. The most used 
are: 

The Free Energy : for processes occurring at constant pres- 
sure. 

This is defined as F = H-TS, and was called the thermo- 
dynamic potential by Gibbs. 

The Available Energy : for processes occurring at constant 
volume. 

This is defined a sA = B-TS and has occasionally been 

termed the free energy by older writers on thermo- 
dynamics. 

Planck s Thermodynamic Potential <j>. 

This is simply defined as <f> = FIT and has certain formal 
advantages in some thermodynamic equations. 

Owing to the fact that these functions have been differently 
named by various authors, it is convenient to remember 
that the quantity F measures the total energy intake H 
minus the heat T?\S absorbed in a reversible process occurring 
at constant pressure. The net energy ‘freed’ in the process 
and available for mechanical work is thus ~F. Similarly, 
the quantity -A is the net energy ‘available’ for mechanical 
work m a reversible process at constant volume. 



II 

THERMODYNAMIC RELATIONSHIPS 


Before discussing the problem of physicochemical equili- 
brium it is useful to consider the minimum number of 
variables such as pressure, temperature, concentrations of 
the various components, etc., required to define completely 
any system. Qualitatively the problem is solved by the 
Phase Rule, symbolized as F — C— P+2, where O is the 
minimum number of molecular components required to 
describe the composition of every phase, P is the number of 
phases, and F is the number of degrees of freedom, i.e. the 
number of variables which can be altered independently with- 
out causing disappearance of one of the phases. 

The quantitative treatment depends on certain relation- 
ships between the thermodynamic functions, and other 
variables such as temperature and pressure. It is particularly 
important to summarize such relationships as are used for 
calculating the values of thermodynamic functions from 
quantities which are more accessible to experiment. 

The heat content and the internal energy. 

For these two functions, by definition 

t x 

= G P H T -H Ti = j Cp dT,' 

T 2 



where G P and C v are the molar specific heats at constant 
pressure and constant volume respectively. 

The magnitude of H T and F Ti relative to suitable refer- 
ence states T 2 can thus be determined by integrating the plot 
of specific heat against temperature, usually graphically. 
Since extrapolation to absolute zero (cf. p. 13) is quite 
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reliable for C v , it is possible in many cases to use the heat 
contents at absolute zero as standards of reference for the 
individual substances taking part in a reaction. This has the 
further advantage that the heat of reaction at absolute zero 
can be directly related to linking energies of chemical bonds. 

Since specific-heat measurements are primarily carried out 
at constant pressure, the values of H are much easier to 
evaluate than values of E , particularly for solids and liquids. 
The values of C v and E are, however, of particular im- 
portance for the theory of the solid state (cf. p. 61). Various 
transformation formulae between the two functions have to 
be used in calculating E from H. These can be derived as 
follows : 


At constant pressure H = E~\-PV and 


dT = 


dT+PdV , 


1 

Also, since E is a function of T and F, 

'0JE\ ( 8V\ 

4 t\8T) x 


dV, 


B. 


- ° r+ [w) T (^) P - 

Combining formulae (1) and (2): 


(i) 


( 2 ) 


For purposes of calculation this expression still has the 
disadvantage that ( dE/dV) T is not usually available, unless 
the equation of state is known experimentally. For a perfect 
gas this term is of course zero. A further transformation can 
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be obtained by using the entropy function, i.e. by intro- 
ducing the second law of thermodynamics. 

Thus in any reversible change cLE = TdS — PdV, 

i.e. m= T m- P 




8VJ T 
W\ (d& 

8tL\8V 


(see equation (1), p. 14), 


so that 


Cr>-C v 


/P 1^7 V 


Still another transformation follows from the fact that P 
is a function of V and T so that 

dP -(§l dT+ Q/ r ' 


or at constant pressure 


(£E 


\8T) v '\8V) t \ 8T) p - 

Eliminating ( BPjdT) v from (1) and (2), since it is a quantity 
difficult to determine experimentally for condensed phases 
(solids and liquids), 

n-c^- T m 


* y \dV ) T \oT 

or since the compressibility 

K = — 

V \dP) T 

and the cubical coefficient of expansion 


P V\dT]p 


r r -^ v t 


This formula is occasionally written with the linear coefficient 
of expansion a = /}/3, in which case 

n-Cr=^-T. 



12 


THERMODYNAMIC RELATIONSHIPS 


The object of these formulae is in every case to calculate 
important quantities from results which are fairly readily 
obtained experimentally. When even the compressibility is 
not known, e.g. at low temperatures, it is possible to use 
empirical formulae such as 

C P -C r = ATi or g£=l+^: 

where the constant A has in each case to be evaluated at a 
temperature where C P and C v are both known. 

Other relationships between H or E and the different 
variables available may be readily derived (cf. p. 16) but are 
seldom used experimentally. 


The entropy S 

The exact evaluation of the entropies of chemical com- 
pounds has been greatly developed in recent years, owing to 
its importance in the calculation of equilibria. 

By definition, if dQ is the heat absorbed in a reversible 
change, the increase in entropy of the system 



dS 

Thus 

(dS) P 

and 

(St-St)p = J - 

Similarly, 

(S Ti ■ 


dQ 
T * 

C P dT 


T 


T x 


T x 

■/< 

T s 


For the calculation of entropies the reference st&te is 
usually chosen to be the crystalline solid at absolute zero, 
since this choice has important theoretical advantages (cf. 
p. 27). From the practical standpoint, however, it should be 
noted that many solids have specific heat anomalies at low 
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temperatures. The contribution of such anomalies to the 
heat content, in a graphical evaluation of the integral 

j G P dT = H, 

is insignificant, of the order of a few calories, but their con- 
tribution to the entropy cannot be neglected. In the plot of 
C P against In T, which is used in a graphical evaluation of 
the entropy, it is not the magnitude of dQ but of dQjT which 
is relevant, so that quite small heat changes make important 
contributions to the entropy when the temperature at which 
they occur is sufficiently low. These specific-heat anomalies 
are discussed in more detail below. 

Important relationships between entropy and other vari- 
ables may be summarized: 

(a) Dependence, on volume and pressure. In an infinitesimal 
reversible change the increase in internal energy dE is equal 
to the heat absorbed minus the external work done by the 
system ; i.e. 

dE= TdS-PdV, 



Similarly, using the relation 

dH = dE+PdV+VdP 
= TdS+VdP , 



Differentiating (A) with respect to T at constant volume, 
(3P\ _m ^ (dS\ d*E 
\dT) v dTdV + \dV) T dTdV’ 
and since ^ Cy 

[dTJr T \8T) v T’ 

8 Z S 1 (8Cr\ 1 8 2 E 

8T8V T\8V ) T T 8T8V’ 

so that the second and fourth terms cancel. 
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Ttos (§i - <*> 

and equation (A) transforms to the 'thermodynamic 5 equa- 
tion of state of a substance: 


Similarly, using the heat-content equation (B), 

(as\ fan 


\8PJt WJr 

and the second 'thermodynamic 5 equation of state is 

xt mtm , fern 


These theoretically exact equations of state are useful 
checks on empirical equations, e.g. in calculating the Joule- 
Thomson effect ( 8H/dP) T in terms of the equation of state. 
They are inserted here since they also permit thermodynamic 
functions to be evaluated when the equation of state is 
known. For a perfect gas, for example, the equation of 
state is TV = RT , 


| for a perfect monatomic gas. 


dP + l-^L 

\8T, 


= ~§dP 


5RdT 
* 2 T 3 


jRln— ~f-ln AL, 


where A is a constant (the entropy constant of a perfect 
monatomic gas) whose evaluation requires statistical 
mechanics, and is discussed later (cf. p. 96). For a perfect 
diatomic gas the entropy due to the translational specific 
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heat has the same value, but terms have to be added for the 
vibrational and rotational specific heats, since these also 
occur in the integral 

(S' = J CpdkiT. 

The problem is further discussed on p. 100. 


The free energy and the available energy 

These functions are defined by the equations 
F = H—TS 
and A = E—TS 


and are more convenient than the entropy in dealing with 
problems of chemical equilibrium. 

Differentiating, 


and 



Similarly, 


F—H 

T~' 



(using equation (A'), p. 14) 


and 



A—E 
T ' 


( 1 ) 

( 2 ) 

(3) 

W 

(5) 

(6) 


Equations (4) and (6) are known as the Gibbs-Helmholtz 
equations, and are of great importance experimentally in 
calculating the variations of F and A with temperature. 

Equations (2) and (5) may be used for calculating the 
values of F and A from equations of state, e.g. for an 
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imperfect gas. 
PV = RT , 


For a perfect gas, whose equation of state is 


'RT 

(dF) T = dP, 


F x -F 2 = RT In PJP 2 
= RTlnCJOt, 

where C\ and C 2 are concentrations expressed in suitable 
units. Similarly, 

A 2 -A 2 = RThiV 2 IV 1 = RT\nPJP 2 . 

Both these changes refer to processes carried out at constant 
temperature. 


Thermodynamic functions and other variables 


Since each of the thermodynamic functions refers only to 
the actual state of the system, and not to its past history, 
each is a complete differential. Mathematically this implies 
that a whole set of analytical relations can be obtained by 
using any pair of the quantities P , F, T , S , 1 E, etc., as 
independent variables. Only the more important relations 
used in calculating the functions from experimentally ac- 
cessible data have been given above, but the general principle 
may be illustrated for the quantities E, S, F. 

Taking any two independent variables x and y, since 
dE = TdS — PdV, 


dx 8x dx 5 


( 1 ) 


dy By By ' 

Differentiating (1) with respect to y. 

By dx By Bx By Bx By Bx By~Bx 
and (2) with respect to x , 

B BE BT BS . m B BS BP BV ^B BV 


(2) 

(3) 
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The property of perfect differentials to be used is that the 
order of differentiation does not matter, i.e. 

d_dE __ 8_dE m 
By Bx ox By * 

whence from (3) and (4) 

dT 38 3P3V 3T 38 3P 3V 

— — » — .. ■— — — m \ O I 

3y cx By ox 3x By 3x By 

Any pair of S , F, T , P, etc., may now be chosen as the two 
independent variables x and y , for which, furthermore, since 
they are independent of one another, dxjdy -may arbitrarily 
be put equal to zero. 

Thus when x = 8, y — F, 

BT\ _ __ (BP 
8V] S ~ \3£ y 

the change at constant entropy referring to adiabatic pro- 
cesses. 

When x = 8, y = P, 



When x = T, y = P, 

3V\ _ _(d8\ 

dTj P ~ {dP) T - 

When x — T, y = V, 

8S\ _ (8P\ 

8V) t ~ \8T) V 

This particular set of equations is known as Maxwell’s four 
thermodynamic equations, but by choosing other initial 
equations and going through the same procedure, a large 
number of relationships between the partial derivatives of 
these functions can be derived. Many of these are mathe- 
matical curiosities with no practical application. 

4360 
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Thermodynamic functions and concentration of the 
components of a system 

The dependence of the functions S, A, F, etc., on the con- 
centrations of the various components which make up a 
system is important for many calculations. It is convenient 
to define Partial Molal Volume, Partial Molal Heat Content, 
etc., of the system by the equations 


dV 

8n x 




^ = H X ~ 2 = H 2 ..., 

8n x dn 2 

where n v n 2 , etc., are the number of moles of each component 
in the system. In general, since V is an extensive function, 


dV = + V 2 dn 2 ..., 

with similar equations for the other functions, so that the 
result can always be written V = n^, H = n x H v and 

all the equations derived for the thermodynamic functions 
of the system as a whole can be split up into formally similar 
equations for the partial molal quantities, e.g. 


F x = H 1 -TS 1 . 

When these partial molal equations are multiplied by the 
number of molecules of the component in question, and 
added, they give the corresponding equation for the system 
considered as a whole. 



Ill 

THERMODYNAMIC FUNCTIONS AND PHYSICO- 
CHEMICAL EQUILIBRIUM 

In principle each of the functions S, A, F, <f> can he used as 
a criterion of thermodynamic equilibrium, though the most 
suitable function to use depends on the physical restrictions 
on the system. The underlying idea is that in any system 
small local fluctuations in temperature, density, and energy 
are constantly arising, owing to the motions of the mole- 
cules, but that when the system is in equilibrium these 
fluctuations do not lead to a spontaneous change in the state 
of the system considered as a whole. The only direction in 
which a spontaneous change could occur is at once given 
from the second law of thermodynamics, 

If the entropy is considered, this can only increase or 
remain unchanged, according to whether the change con- 
sidered is irreversible or reversible. Formally this may be 
stated as follows: 'Keeping the volume and internal energy 
of the system constant, small fluctuations in the other 
variables must lead to a zero change in entropy if the system 
is in equilibrium, and to a positive change in entropy if the 
system is not in equilibrium 5 , i.e. 

SS ^ 0. 

Similarly, at constant heat content small fluctuations must 
be subject to the condition SF < 0, or at constant internal 
energy to the condition S A ^ 0. If it is possible to express 
the dependence of the thermodynamic functions of a system 
on the concentrations of the components, etc., then the con- 
centrations at equilibrium can be obtained by equating the 
first differential of the appropriate function to zero. 

From a purely formal standpoint, the choice of which 
independent variables to use in expressing the thermo- 
dynamic properties of a system is summarized in a rule 
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given by Massieu. This states that the thermodynamic func- 
tion whose maximum or minimum, at constant value of the 
independent variables, expresses the condition of equilibrium, 
gives the thermodynamic properties of the system by simple 
differentiation, i.e. in the most elegant form. The rule is 
best illustrated by a table. In the third set the results are 
described for the function <f>, since those for the function F 
are closely related and have already been described. 

Independent variables. Equilibrium conditions. Characteristic function . 


(These have to be kept 
constant in calculating 
the effect of fluctuations 

in the dependent 
ables.) 

vari- 

E and V 

SS = 0 

(i.e. no heat exchange with the surroundings) 

giving (cf. p. 17) 

(BS\ 1 (dS\ P 

\dE/r ~~ T \dV/ e T 

V and T 

8A < 0 

giving 

(dA\ _ (8A\ __ 

XdVjT \8tJv 

P and T 

8<f> < 0 

(E+PV) 
~ S+ T 

(H\ _v h 

\dP ) t T \dTJp T 2 * 


In connexion with this function it may be noted that 
E -~ T { T ^ +P %\ 

-S = t + T S ±, 


so that the thermodynamic behaviour is completely determined when <f> is 
known as a function of T and P. 

For practical purposes the most important function is the 
molal free energy, or for systems involving more than one 
component, the partial molal free energy. Before discussing 
how this is used for calculating equilibrium constants the 
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proof may be given that in a system at equilibrium the free 
energy /unit mass in each phase must be the same. 

The condition for the coexistence of two phases of a sub- 
stance in equilibrium is SS — 0, when the subordinate con- 


ditions are 


Sm = 0, . 8E 


SV = 0. 


Suppose the fraction x to be in phase 1 , so that if m is the total 
mass of the system, the mass in phase 1 is mx and in phase 2 
m( 1— #). If S ± and S 2 are the entropies per unit mass in 
the two phases, the total entropy of the system is S, where 
S/m = #^+(1 — x)S 2 . 

The total volume is F, where 


F/m = x V x +(l — x)V 2 . 

The total energy is E , where 

E/m = xE x +(l — x)E 2 . 

Applying the condition SS = 0, since m is constant, 
x SS-l + (1 — x) SS 2 -(- ( S 1 — S 2 ) 8x = 0. 


Now 


8S X = 


T x ’ 


S 8 , = 


S E z +P 2 8V 2 

T 0 


Eliminating SV 2 and 8E 2 by means of these equations, 




Since the variations SE ± and 8V ± are independent, this equa- 
tion can only be generally true if the intensive functions are 
the same, viz. T x = T 2 , Pi = p» and furthermore if 


-S x . 


Ei+PV x Q , E 2 +PV 2 


T 


$2 + 


T 


i.e. F ± = F 2 . 

In the case of systems involving more than one component 
a similar proof shows that the partial free energy/unit mass 
must be the same in every phase. This quantity has therefore 
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been termed the thermodynamic potential of the component 
in question. If the molecular weight is chosen the same in 
every phase, the partial molal free energy must also be the 
same in every phase. 

Finally, it must be noted that the definition of thermo- 
dynamic equilibrium rests upon the criterion of the effect of 
small fluctuations on the thermodynamic functions. This 
leads to a physically satisfactory treatment of equilibria 
such as the vapour pressure of solid nitric oxide, or the 
equilibrium 2H 2 +0 2 ^ 2H a O, 

though if sufficiently large energy changes could take place 
within reasonable time the nitric oxide would change com- 
pletely into nitrogen and oxygen at low temperatures, and 
the hydrogen and oxygen nuclei would probably change into 
helium atoms. Practically all thermodynamic equilibria are 
pseudo-equihbria in the sense that it is not the absolute 
minimum of free energy with respect to all possible changes 
which is discussed, but only with respect to changes which 
could take place within reasonable time. 

Standard free -energy calculations and the equilibrium 

constant 

For practical purposes the most useful criterion of equili- 
brium is BF = 0 , since most experiments refer to processes 
taking place at constant pressure. Its application may be 
illustrated by calculating the free -energy change in a reaction 
between gases, such as C0+H 2 0 ^ C0 2 +H 2 . Neglecting 
the correction for departures from the gas laws, which is 
irrelevant for the present purpose, the tree-energy change 
on changing the concentration of any one component is 

F x ~F 2 = XTlaOJOt (cf. p. 16). 

This expression may be formally simplified by choosing for 
the standard state of reference the free energy F 2 at unit con- 
centration (or one atmosphere pressure where partial pressures 
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are used), i.e. C 2 = 1. Writing the free energy in this standard 
state as F 0 , the free energy in any other state is given by 

F x —F q = BTlnC ly 

and the total free -energy change in the reaction is 

AF = 'ZF 1 -ZF 0 = 

So far the C ± ’s have been left arbitrary; by selecting them 
to be any set of equilibrium concentrations for the mixture, 

2-^1 = 0 and FT^laC e(l = — 2-^0- 
This follows from the fact that for an infinitesimal transforma- 
tion of 8n molecules in the equilibrium mixture 

S^(C0+H 2 0) -> Sn(C0 2 +H 2 ) 

the tree-energy change 

8n 2 == $F = 0 

at equilibrium, so that 

2-^leq = 

Thus for any set of equilibrium concentrations 
-2^0 = BT^laC Q(l = RT]nK } 
and since the left-hand side is a constant, 

K = constant, 

where the equilibrium constant K is called the mass-action 
constant. The magnitude of the mass-action constant in terms 
of the standard free energies of the individual components is 
thus at once obtained in this method. 

This equation is of fundamental importance in the applica- 
tion of thermodynamic functions to the calculation of equili- 
brium constants. It would, however, lead to no great prac- 
tical progress if the only way of determining the thermo- 
dynamic functions experimentally were from studies on 
equilibrium constants, the electromotive force of reversible 
cells, and the limited number of other physicochemical 
measurements from which free energies can be calculated. 

A much more general means of calculating thermodynamic 
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functions is from heat data, which can be obtained fairly 
easily in a large number of cases by direct calorimetry. The 
Gibbs-Helmholtz equation 

(3F\ q F-H 

~T~ 

gives on transforming and integrating 

J = ^ = +i, 

where I is an integration constant. (The same result is 
obtained even more directly by integrating the equation for 
cf> on p. 20.) Thus values of F can be calculated from values 
of H provided the integration constant I is known. Since 
F = — RThxK, 

an alternative way of writing the above equation directly 
in terms of the equilibrium constant is 

^ K -\^ dT - 1 ■ 

Using this formula, interpolation formulae for the values of 
F or K at different temperatures can be constructed from 
the following experimental data: 

1. A knowledge of the value of H over the range of 
temperatures for which the interpolation formula is to be 
constructed. Since in any reaction (dH/3T) B = Cp (cf. 
p. 9), it is always possible to write 

^ = ^+{( 2 ^) dT, 

T x 

so that it is sufficient to determine experimentally the value 
of H at any one temperature, and the specific heats over 
the range of temperatures for which the formula is required. 

2. In order to determine the integration constant /, the 
value of F at any one temperature is also required, and must 
be measured experimentally. 

An alternative way of stating the problem of calculating 
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chemical equilibria from heat data may also be given, since 
it is useful in what follows. In any chemical equilibrium 

F = —RThiK = J n>i£, 

where the suffix denoting that the free energies refer to the 
standard value at unit concentration is omitted, since no 
confusion arises. 

This may be transformed into 

F = % n r (H r -TS r ) 

= 2 n r H r — T 2 n r S r 
= H-T^n r S ri 

where H is the change in heat content of the system. Thus 
the evaluation of F again resolves into an evaluation of H and 
of the individual entropies 8 r of the substances taking part 
in the reaction. Since any entropy may be written 

Tr 

S Ti = S Tt + J Op d In T, 

T s 

the construction of an interpolation formula for J 7 as a 
function of temperature is again seen to require 

(1) The determination of II at any one temperature. 

(2) A knowledge of the values of C P over the whole range 
of temperatures required. 

(3) A determination of the individual entropies at any 
one temperature. 

In a very large number of cases the experimental diffi- 
culties in evaluating the thermodynamic functions arise in 
determining at least one value of the free energy or entropy 
change in the reaction, since some sort of a reversible process 
has to be constructed and studied. It was with a view to 
eliminating the necessity for this measurement that the third 
law of thermodynamics was developed. 
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The evaluation of thermodynamic functions from heat 
data 

The practical calculation of the functions H, S 9 F s etc., 
from calorimetric measurements alone rests upon two hypo- 
theses first enunciated by Nernst,f and sufficiently verified 
and clarified since to be incorporated as the Third Law of 
Thermodynamics . 

The evaluation of H: Nernsfs first hypothesis. 

The calor im etric determination of the change in heat con- 
tent for a chemical reaction does not offer major experi- 
mental difficulties, provided a suitable temperature can be 
chosen for the measurement. The value of H at any other 
temperature can then be calculated from the equation 

Bt-Bt, = f (2 0 P ) AT. 

T s 

The value of H at absolute zero is of particular interest. In 
order to obtain it, it is necessary to extrapolate specific heat 
measurements below the lowest temperature of measurement. 
Nernst assumed that for 'condensed 5 phases (solids and 
liquids) 

lim 

which means that the extrapolation from such low tempera- 
tures as can be readily obtained, to absolute zero, gives quite 
reliable values of the integral 

T 

J Cp AT. 

0 

This is readily seen from Fig. 2. The theoretical expression 
C v = ocT z obtained by Debye (cf. p. 66) is often used for the 
low-temperature extrapolation. Keesom has shown that 
for some metals C v = ocT is closer to the experimental 

t Gottinger Nachrichten (1906), 1. 


— lim Cp 


T—>0 
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results, f but the exact expression used makes little difference 
to the area under the curve, since (dHjdT)jp is very small in 
either case. 



Pig. 2. 

Evaluation of 8; use of Planch's formulation of the heat 
theorem. 

The evaluation of absolute entropies (i.e. entropies re- 
ferred to the crystalline solid at absolute zero as standard) 
is a much more delicate problem than the evaluation of 
heat contents, and is one of the principal aims of contem- 
porary thermodynamics. From 8 and H the value of F 
and of the equilibrium constant is readily computed. By 
definition, 

Sp = ^CpdlaT +S 0 , 

where S 0 is the entropy at absolute zero. 

Planck first suggested that the entropy of crystalline 
bodies at absolute zero was probably zero, i.e. 8 0 — 0. A less 
sweeping statement is that every body has a finite entropy 
which is positive but which may become zero at absolute 
zero for perfect crystalline substances. This form of stating 

f Physikalische Zeitschrift , 35, 939 (1934). 
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the heat theorem has the great advantage that it is much 
easier to see its statistical significance. Its connexion with 
'Nemst’s two original hypotheses 


v 8H A 
hm-= = 0 

T=0 81 


dF 

and lim— = = 0 

t=o8T 


is readily seen, since 

8F _ 8H_ t 8S_ _ m_ a_Cp 
8T~8T 8T~8T T ’ 

which only tends to zero at absolute zero if 
S 0 = 0 and lim G P — 0. 

The equations are written in the form of a limit to indicate 
that the approach to zero is asymptotic, i.e. that C P decreases 
with T with a power of T greater than unity. 

When the heat theorem applies, the calculation of entropies 
from experimental data merely involves an evaluation of the 
integral 

$ = J C P dlnT, 

since 8 0 = 0. The difficulty of specific-heat anomalies at 
very low temperatures, in upsetting experimental evalua- 
tions of entropy, has already been referred to. The most 
direct test of how far such anomalies are widespread is to 
see how far the heat theorem applies to practical cases. It 
is useful to note that whenever an anomaly occurs at so 
low a temperature as to escape measurement, the entropy 
evaluated from specific-heat measurements will be too low, 
and can only be made to agree with the entropy evaluated 
by other means by giving S 0 a finite positive value. A number 
of apparent departures from the heat theorem can be 
explained in this way. 



IV 

APPLICATIONS AND VERIFICATl^fe^fc*^ 
HEAT THEOREM 

The method of applying the heat theorem in practice is best 
grasped from the methods used to test its validity. 

(a) V erifications for condensed systems . 

Whenever the equilibrium constant in a system involving 
only solid or liquid phases can be measured directly, a means 
is available for checking the heat theorem. The number of 
cases is not very large. One example is the calculation of 
the transition or equilibrium temperature between two solid 
phases. The method is simply to determine experimentally 
one value of H and the heat capacities over the whole range 
of temperatures. The absolute entropies are then calculated 
assuming Planck’s law, and from the heat content and 
entropy the free energy and equilibrium constant are readily 
found. The result should agree within the limits of experi- 
mental error with the value determined directly. 

Thus grey tin is an allotropic modification of white tin, 
and is less dense. Below 18° C. white tin tends to change 
spontaneously into grey tin, and above 18° C. grey tin tends 
to change spontaneously into white tin, this being the transi- 
tion temperature under ordinary pressures. Prom specific- 
heat measurements extended to very low temperatures, 
assuming S 0 = 0 for both forms, 

$ 25 °c. = 11*17 (white tin), $ 2 5°c. — 9*23 (grey tin ). 
Difference + 1 * ! 94. 

(The entropy of the less stable form is necessarily smaller, 
since a spontaneous change at constant volume can only take 
place with increase of entropy.) 

This result has to be compared with the directly deter- 
mined value of A/S^ ey ^ white . A H is measured calorimetrically. 
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e.g. by dissolving each allotrope in acid in a calorimeter. 
A F is measured from the e.m.f. of a cell (£, 

grey tin | electrolyte | white tin, 

since the free-energy change for the passage of one faraday 
of electricity F is 

— N(£F = AF = —RTlnK. 

From the equation A F = AH—TAS the experimental 
value of AS = +1*87 (grey -> white) is obtained. This com- 
pares very favourably with the value directly obtained from 
heat data. 

Other systems involving only reactions between solids, 
which have been investigated, are 

sulphur (rhombic) ^ sulphur (monoclinic) 

Pb+I 2 ^ Pbl 2 
Pb+2AgCl PbCl 2 +2Ag 
Hg+AgCl^HgCl+Ag. 

The agreement in each case is within the limits of experi- 
mental error. 

(b) Verification from condensed systems in equilibrium ivith 
the gaseous phase. 

The heat theorem applies directly only to ‘condensed’ 
systems, and though at very low temperatures and high 
densities even gases probably show degeneracy, with C B 0, 
no experimental measurements have been carried sufficiently 
far to make use of this fact in the evaluation of S. In conse- 
quence it is not possible to evaluate the absolute entropy 
of a gas using only thermal data and the laws of thermo- 
dynamics. 

Statistical calculation (cf. p. 98) shows that the entropy 
constant of a monatomic perfect gas in the expression 
obtained on p. 14, 

5 = 
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is given by 


^ __ { < Z'Trmk)*lc 


where m is the mass of the molecule, k is Boltzmann’s con- 
stant (k = RjN), and h is Planck’s constant. 

In terms of the atomic weight M this expression transforms 
to S = R\a(T%M%)jP — 2*300 cals./° per mole, where the unit 
of pressure is the atmosphere. This statistical result may be 
used in a number of ways to verify the theorem that S 0 = 0 
for crystalline substances. 


I. The vapour -pressure equation for solids 

The most direct way of verifying the heat theorem is from 
the vapour-pressure equation, either for solidified inert gases, 
or for metals such as Hg, Cd, Zn which give monatomic 
vapours. 

On vaporization the entropy increase per mole is 

as = s gas -s 80lid = A 

= [Pln(?W*)-.KlnP— 2-30]- J j C Pi din P+P 0 ], 

where A is the latent heat of vaporization per mole, C Pi the 
specific heat of the solid, and S 0 the entropy of the solid 
at 0°K. 

This expression is merely another way of writing the 
familiar vapour-pressure formula of a solid giving a monato- 
mic vapour, since it transforms to 

hlP= J Cp ^ lnP+ [ fln - M -^“§]’ 

. . , 5RT 

or smce A — A 0 -] — , 




to 

In P 
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The temperature independent term in this equation is called 
the vapour-pressure constant i. 

In testing the validity of the heat theorem it is usually 
more satisfactory to compare the entropy of the gas at any 
temperature, calculated directly from the equation 

# gas - RkiTmiP, 

with the experimental value determined from the vaporiza- 
tion prOCeSS 0 \rm \ o 

^g-as — Aj J. ~r -Osolid- 

In the experimental determination the entropy of the solid 
is taken as, zero at the absolute zero, in accordance with the 
heat theorem, and the value of A is computed from the slope 
of the vapour-pressure curve. Typical results at 25° C. and 
one atmosphere are 


Gas 

Experimental 

entropy 

Calculated 

entropy 

He 

29-2 

29-8 

Ar 

36-4 

36-7 

Cd 

40-0 

39*8 

Hg 

41-3 

41*5 


The good agreement provides a simultaneous check on statis- 
tical theory and on the heat theorem. 

An older method, avoiding the explicit use of entropies, 
depends on the comparison of the vapour-pressure constant 
with the theoretical value. From the Clausius -Clapeyr on 
equation ^ 

dP (^gas ^solid ) P % 


At the very low pressures to which the experimental data 
refer v gas = RTjP per mole and v soYia < %l S * The equation 
transforms to 1 > 


P dT ~~ RT 2 ’ 


j 


A dT , . 

RT 2 + % 


or 


In P 
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where i is the vapour-pressure constant. In this expression 

(Cpfeas) ” ^P(solid) ) 5 

i.e. A = A 0 + J ( C g — C s ) dT , , 

so that lnP = - A+ JJ *L(C g -C s )dT-, 
or, using partial integration, 

lnP= -^p +5 J ^- c ^ dhxT -^rj ( c„-c.)dT+i . 

This may be compared with the direct expression for the 
entropy of vaporization: 

s eas -s aom = A = ^+1 J (C7,-qy <*p, 

Sgas = J CJ, din T-B In P+31, 


^solid == J Cg C^ln T + $ 0 . 

Rewriting the equation for the entropy of vaporization, 

+ J (0,-CJdlnT-M, 

from which it follows that provided the same units of 
pressure are used, and provided the heat theorem applies 
for the solid phase (S Q = 0), 



(Some authors include the term — ^ J C g dT in the value of i, 

since for a perfect monatomic gas it is independent of 
temperature and equal to |.) The theoretical value of the 
vapour-pressure constant is simply calculated from that of 
the entropy constant for the appropriate vapour, and is 
compared with the experimentally determined value. 
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II. Comparison of equilibrium constants with vapour - 
pressure constants 

Although the use of free energy and entropy functions 
provides the most compact and general method of calculat- 
ing equilibria, the vapour-pressure constants or 'theoretical 
chemical constants 5 have been frequently used in the past, 
particularly by German authors. The relation between the 
two different methods of calculation may be illustrated for 
the reaction aA+bB+cC ^ dD+eE.l. . 

A F = — RTlnK = AH—TAS 

T T 

= A H 0 + J (2C P )dT-Tf (ZC^dlnT-TZS, 


0 (solids)" 


0 

T 


^ J ( 2 Qpteas)) dhlT T 2 $0 (gases) • 
o 

This expression may be rewritten 
-In# — AH ° 




C P dT 


RT ' RT 

i f 2 (Csoiids+C'* ) din T— y 


R 


O(solids) 


2 ^0 (gases) 

R 5 


j R 

where S 0 refers to the entropy at 0° K., H 0 to the heat content 
at 0° EL. 

This may be compared with the reaction isochore, 
dln.fi: AH 


dT " 

which gives on integration 


" RT 2 ’ 
RT 2 + 


I. 


Substituting for AH, 

In K^- AH ° 


RT ' R 

or integrating by parts, 

AH, 

RT 


]nK 


1° L f y 

V RT J Z 


IfijsfCZOJdT+I, 


Op dT + 


dT 

T 


+ !■ 
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From a comparison with the previous expression it is seen 
that the integration constant 


r V solids) 

Z, R 


+ 2 


S, 


'OCgases) 

R ’ 


so that if the heat theorem applies for solids , making 

2 ^O(solids) == 

the integration constant can be expressed in terms of the 
individual vapour-pressure constants, according to the 
equation „ « 

j __ ^ v i = 2 ° 0(gases) 


(v is the number of molecules of the gas participating in the 
reaction). 

It is thus possible to test whether the heat theorem applies 
to solids by comparing the calculated values of the integra- 
tion constant with the sum of the individual vapour-pressure 
constants. Although heterogeneous equilibria are not so 
suitable for exact calculation as the vapour pressure of a 
single substance, on account of experimental errors, this 
method has the theoretical advantage that no appeal is made 
to statistical theory, since no use is made of the theoretical 
entropy constants of gases, provided the experimental 
vapour-pressure curves are known. This method is therefore 
directly applicable to equilibria involving polyatomic gases. 
Typical results may be listed as follows :f 


Reaction 

j Limits of I 

j {from reaction 
j equilibrium ) 

Limits of 2 
(from vapour 
pressures) 

3H 2 +N 2 ^ 2NH a 

-7-04±0-10 I 

-8-34±0-09 

N 2 +Oo ^ 2NO 

0-93±0-28 

0-61±0-15 

2C0 + 0 2 ^ 2CO a 

— 0-8±0-25 

-l*38±0-20 

H 2 +I 2 ^ 2HI 

-1-5±0*12 

-2-36±0*22 

2CO ^ C+C0 2 

— 0-86±0-18 

-1-01±0-18 

-Hg (vap ) + 0 2 ^ 2HgO 

4-32±0-18 

4-44±0-09 

CaO + C0 2 ^ CaC0 3 

0'9±0*15 

0-91±0-15 

H 2 +HgO ^ Hg+H 2 0 

-3-63 + 0-16 

-3-69±0-03 


f Eucken and Fried, Zs. f. Phys. 29, 36 (1924). 
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In some cases the agreement is seen to be within the limits 
of -uncertainty, which are rather large in this method of 
verification. Reactions involving hydrogen nearly always 
give very poor results. For reactions involving CO and NO 
the uncertainties of this method are too great, but a more 
precise method (cf. p. 74) shows that the discrepancy is real, 
and gives some clue as to its origin. 

Theoretical chemical constants of polyatomic gases 

Owing to their role in the integration of the reaction 
isochore, the vapour-pressure constants of the various com- 
ponents have been given the name of theoretical chemical 
constants. Once their value is known it is in principle 
possible to calculate all chemical equilibria from heat data 
alone. 

The theoretical chemical constants of gases and the closely 
related entropy constants can be calculated from statistical 
theory. For monatomic gases the entropy constant 21 in the 

expression S P = B ln^- 21 is 21 = ^ 7rm ^ ^ e ~ . As is discussed 

later (p. 54) this entropy constant has to be multiplied by 
various factors if the atom or molecule of the gas is capable 
of taking up energy in other forms than translational energy, 
provided such activation is appreciable over the region of 
temperatures where the equilibrium is discussed. For the 
present purpose these rules may be summarized as follows: 

(1) If the atoms have an electronic multiplicity j > 0, 
21 has to be multiplied by the factor g = 2j+l, since in the 
modified statistical calculation 2j-\~ 1 atoms with closely 
similar but still distinguishable energy levels can now exist 
in the gas phase, for every atom which is there when the 
multiplicity is unity. This has been experimentally verified 
in the case of thallium, which has a ground state 2 P^ and 
g = 2, and a correspondingly larger entropy of vaporization. f 
t Coleman and Egerton, Phil . Trans. 234, 177 (1935). 
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(2) For polyatomic molecules the term has to he further 
multiplied by the factor for rotational energy 

87 t 2 IJc 1 

The moment of inertia I — %J(I ± J 2 / 3 ) of the molecule has 
to be evaluated from the separation of rotational bands in 
the band spectrum. The symmetry factor s only applies for 
symmetrical molecules such as N 2 , 0 2 for which it is 2 
(cf. p. 120). 

Practical calculation of the position of equilibrium in 

reactions 

Before discussing the heat theorem further, the practical 
utilization of chemical constants may be referred to. For 
most purposes departures from the heat theorem do not 
much affect the value of the equilibrium constant calculated 
on the assumption that it holds. The only systematic method 
of calculating equilibria is to determine the entropy and heat- 
content functions for as many substances as possible, since 
from these the values of A F — — RTlnK are at once 
obtained. The use of vapour-pressure constants is formally 
equivalent to this method, but in practice it leads to the 
tabulation of data from which the desired results can only 
be obtained with much more labour. 

As an example of the calculation of equilibria not accessible 
by other means, the case of ionization in a star may be 
quoted. 

In the process A ^ A + +e 

the entropy change may be written 

AS = S A +-\~-8 e — $ a P er mole. 

Apart from changes in multiplicity of the atom on ionizing, 
which affect the entropy constants by affecting the factor g , 
but which will not be discussed here, 

S A + = S A and AS = S e . 

For the equilibria in question the pressures are so low and 
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the temperatures so high that the electron gas is very dilute, 
and may be treated to a first approximation as a perfect 
monatomic gas of atomic weight M = 0*000544, from which 
the entropy constant of the electron gas can be calculated. 
LTsing the formula (p. 31), 

AS = S e = §B]nT+%B1nM—2-30 
= 5 In T— 24*70. 

The change in heat content is approximately equal to the 
ionization energy per mole. If V is the ionization potential 
in volts, 

A H 0 — 23*3 V X 10 3 calories, 

AH = AH 0 + j (C A +C-C A+ ) dT = AH 0 + 5 -ZT, 
since C A — C A +, 

A F = -BT\n~ = 23-3Fxl0 3 -53 r, ln2 1 +29-7T 


, A + e 
log — = 


23-3F X 10 s , 11-5. _ 29-7 

-4-57T +4^7 l0gy -?57 


5*lF X 10 3 

T 


- 2*5 log T — 6*5. 


Since V = 4 to 5 volts for the alkali metals, the right-hand 
side is only positive when T ^ 10 4 °C. In the spectrum of 
the sun the lines for Na are strong, for K weak, and entirely 
absent for Rb and Cs. This is due to the fact that at the 
same temperature the abundance of unionized alkali metal, 
which is required for the spectrum, is largest for the metal 
with the largest ionization potential, as is clear from the 
above formula. 


Approximation formula due to Nernst 

For very rough calculations of the position of equilibrium 
in a reaction an approximation formula due to Nernst often 
gives quite reasonable estimates. This may be written 
A/7 

logZp = - + 2 »(i-75iogr)+ 2 C, 
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where products are as usual accounted positive in the summa- 
tions, and n refers to the number of molecules of each species 
taking part in the reaction. Tor reliable results 2 > ± 1. 
C is an empirical ‘chemical constant’ quite distinct from 
the vapour-pressure constant and only for use with this 
empirical formula. The Moving list has been given by 
Nemst.f 


h 5 

1-6 

HC1 3-0 

H 2 S 3-0 

n 2 

2-6 

HBr 3-2 

SO, 3-3 

o 2 

2-8 

HI 3-4 

CO, 3-2 

Cl 2 

3-1 

NO 3-5 

NH S 3-3 

Br 2 

3-2 

CO 3-5 

H,0 3-6 

I. 

3-9 

N,0 3-3 


For solids C 

= 0 and for atoms 0 = half the value for the 

corresponding molecules. 

An example of 

how far equilibria 

can be represented by this formula is the 

reaction 

H 2 +C1 2 ^ 2HC1, 

A#i5°c. — ~ 

-44,000 cals. 



log K P 


Temperature 

obs. 

calc. 


300° 

33-2 

33-38 


1,000° 

10-43 

10-93 


2,000° 

5-77 

6-11 


though the equation should be used with caution, since the 
values of 0 only fit a limited set of data. 


f The New Heat Theorem, 
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According to the Nernst Heat Theorem as formulated by 
Planck, the absolute entropy of crystalline solids S 0 — 0 
when T = 0. As formulated by Nernst, 


v 8H A r 8F 
i. im ^ = 0 = hm — 


for condensed systems. 

The chief thermodynamic methods of verifying this 
theorem have already been referred to, and involve equilibria 
with one or more components. In principle, however, any 
method of studying the change of these functions with 
temperature, as the absolute zero is approached, will serve 
to test the theorem. An outstanding example is liquid 
helium. 


Liquid helium 

The properties of liquid helium are of considerable thermo- 
dynamic interest, since it is the only substance which remains 
liquid down to the lowest temperature attainable, and its 
thermal behaviour shows that the heat theorem is obeyed. 

Surface tension. The change in free energy associated with 
surface changes is AF = o ACl, where o is the surface tension 
and Q the area, so that the heat theorem would predict 


lim 


dAF 

~w 


0 = lim 


dcr 
dT ’ 


which is what is actually observed, f 

Change in melting-point with pressure . Although helium 
remains liquid down to the lowest available temperatures 
(0*8° K. or less) under its own vapour pressure, i.e. shows no 
triple point, it solidifies under, pressure, the melting-point 
being 4-23°K. at 140*5 atm. and 2-15°K. under 25*3 atm. 
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As T 0 the melting-point curve tends to become parallel 
to the temperature axis. This is what would be expected 
from thermodynamical reasoning. The Clausius-Clapeyron 
equation states that the change of entropy on melting is AS, 
where AS _ dP 

AF ~ dT’ 



and since AF ^ 0 on melting, AS 0 as T -> 0 means that 
dP/dT -> 0 if the heat theorem applies. This is what is 
actually found. 

The fact that helium obeys the heat theorem, in spite of 
its remaining a liquid down to the lowest temperatures at 
ordinary pressures, implies that the liquid has some means of 
decreasing its entropy. This is probably connected with the 
transition into a second liquid form, Helium II, below 2*3° K. 
The transition can be followed from measurements on the 
density, surface tension, specific heat, etc. There is no latent 
heat change when liquid I is transformed into liquid II. 
This could only be expected on thermodynamic grounds if 
two immiscible liquids existed in equilibrium at the transi- 
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tion temperature. There is, however, a sharp maximum in 
the specific heat, similar to that observed in the ammonium- 
chloride effect (p. 75) and other internal transformations 
within a single phase. The slow rise on the low-temperature 
side and sharp fall on the high-temperature side are quite 
characteristic. 

Evidence on the nature of Helium II is obtainable from 
thermodynamic considerations. Since the entropy change on 
melting progressively tends to zero as the temperature falls, 
and the solid has a progressively lower entropy on account 
of its state of order, and progressive decrease in heat content, 
the entropy of Helium II likewise gets progressively lower, 
suggesting it is in a state of quasi- crystalline order. A quali- 
tative description is to compare Helium II with a crystalline 
solid in which the energy of the null-point oscillations ihv of 
the atoms about their positions of equilibrium (cf. p. 61) 
exceeds the work required to shift an atom from one position 
to the next. The 'solid’ thus shows no rigidity to shearing 
stresses, i.e. behaves apparently like a liquid. When the 'solid’ 
is compressed sufficiently, the work of shifting the atoms 
from one position to the next in the lattice increases, and 
a finite rigidity is observed, marking the transition from 
Helium II to solid helium under increased pressure. 

Supercooled liquids (glasses) 

The case of helium, where the liquid form is in true equili- 
brium down to the lowest temperatures, is in sharp contrast 
with the class of supercooled liquids, or glasses, all of which 
fail to obey the heat theorem unless it is properly stated. 

The ease with which a liquid can be cooled below the 
normal freezing-point is greatly increased as the number of 
associative groups such as hydroxyl is increased in the 
molecule. Typical glasses are given by glycerine and ethyl 
alcohol, and by the oxides and acidic salts of the oxides 
^2^5? As 2 0 3 , B 2 0 3j SiO^. 
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Two things favour the passage into a glass: 

(a) The true crystal form is complex, so that the number 
of fresh crystallization nuclei formed in unit time below the 
melting-point is small. 

(b) The rate of growth of the crystal on to a nucleus is 
small on account of the high viscosity. For substances 



possessing these characteristics the properties of the super- 
cooled liquid can be studied with ease well below the melting- 
point. Thus the Heat Content, Volume, Viscosity, Dielectric 
Capacity, and Electrical Resistance for a supercooled liquid 
follow a general curve of the type shown in Fig. 4 . To 
take one property of great importance for glasses, the 
viscosity of the liquid rapidly increases on cooling according 
to a law to 

log t? = — A + - =~— 

up to the temperature at which the glass becomes rigid. Over 

d 2 Pr 

the small interval marked out on the curve v increases 

dT 2 1 

from about 10 3 through the workable range of the glass 
(10 5 to 10 10 ) to the point of rigidity 77 10 13 . It may be 
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noted that a good glass that can be annealed without too 
much difficulty must not have the working range of tempera- 
tures too narrow. Tamman ascribes for example the change 
in the dielectric constant to the gradual rotation of permanent 
dipoles in the glass as the temperature is raised. The change 
from glass to supercooled liquid has thus certain analogies 



with the ammonium- chloride effect (p. 75), except that 
there is* no change in crystalline form, associated with in- 
creasing rotation in the solid. 

The entropy difference between a solid and the glassy 
supercooled liquid can be evaluated by careful specific-heat 
measurements from the melting-point down to very low 
temperatures. Fig. 5 shows the results obtained for glycerine 
by Lewis and Gibson. The initial drop for the solid compared 
with the liquid is due to the loss of heat contributions 
by atoms and groups when these are firmly bound. The 
gradual subsequent drop is due to the Debye oscillations of 
the lattice. The specific heat of the supercooled liquid 
remains approximately 3E/2 above that of the solid, till the 
glass’ is formed at about 180° K.; within a few degrees the 
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specific heat drops to the value for the crystal, since the very 
small differences in the characteristic frequency (cf. p. 64) 
of the crystal and glass are only noticeable at very low 
temperatures. 

From these two curves the quantities liquid and S so i id can 
be calculated at the melting-point. If S 0 = 0 for both solid 
and glass, the equation 

ASeasion = Al 5° n = Sliauld-Ssolid 

should hold. Actually a discrepancy is found in the sense 
that S 0 for the glass must be taken about 4*6 entropy units 
above the value for the crystal, in order to obtain agreement 
at the melting-point. As has been explained on p. 28, such 
a positive value of S 0 is sometimes due to a transition below 
the lowest temperature of measurement making a contribu- 
tion to the entropy. In a glass such a transition is, however, 
very unlikely, and the conclusion is that glassy glycerine 
does not obey the heat theorem. For other glasses which 
have been measured discrepancies of the same order of 
magnitude have been found. This result is not very surpris- 
ing, since from the point of view of thermodynamics the 
transitions crystal-glass or liquid-glass are not really rever- 
sible, and if a molecule escapes from the glass, the structure 
on recapture and energy gained are to some extent indeter- 
minate, since it has not to occupy a definite position in the 
lattice, as it would in a crystalline solid. 

The entropy of mixing: solid solutions and the heat 
theorem 

The general expression for the entropy increase on allowing 
two perfect gases to mix by diffusion, keeping the total 
pressure constant, is readily obtained. For each component 
separately, before mixing, 

&i = n^GjplogT— -KlogPi+fci), 

S 2 = n 2 (C P log T—JtlogP 2 +lc 2 ), 
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where the initial pressures of the %, n 2 moles are Pj = p = j 
and the entropy constants are h v h 2 , For the mixture, 

S = I« 1 (C Pi Iog3’-JJlogPl+y, 
where PJ, etc., are the partial pressures 
p;=%p/(wi+w 2 ) 

i.e. S = % ^logP-PlogQP+y. 

The entropy increase on irreversible mixing by diffusion is 
tks AS = S-(^) = -fi^ilogC, 

Dividing by (%+w 2 ), -FJ CplogCi P er Two 


densed phases are involved. At first sight it would appear 
to be necessary to abandon the heat theorem for mixtures. 
Statistical considerations suggest, however, that the solution 
or mixture may have a number of ways of decreasing its 
entropy with falling temperature, and thus of approximating 
to the heat theorem, though in practice the processes involved 
may occur so slowly that experiments on solid solutions do 
not reveal them. One of these processes is the formation of 
ordered from disordered alloys (cf. p. 79), with a decrease in 
entropy of the solid solution. For gases it should further be 
noted that the mixture formula would have to be modified 


first is known as the Gibbs paradoxon. Gibbs pointed out 


by experiment, there should be no entropy increase on 
mixing, since the case could not then be said to differ from 


are put into communication. In such a case the molecules 
still interdiffuse, spontaneously, but there is no entropy 


lsed for the entropy is then no longer valid. 

Other departures from the heat theorem, such as solid 
lydrogen, carbon monoxide, and nitric oxide, are best dis- 
cussed in terms of statistical theory, and are referred to 
.ater (p. 104 for hydrogen; p. 74 for CO and NO). 


the molecules of the two gases differ for the result still 
to hold? 

As will be seen later, the entropy increase on mixing has 
been confirmed even for mixtures of rotational isomers such 
as ortho- and para-hydrogen; for mixtures of isotopes experi- 
ments are not yet available but probably give the same result. 



A second difficulty arises for solid solutions. Exactly the 
same expression for the entropy increase on mixing holds for 


A# — ^mixture ^ ^components) 

must decrease to zero towards absolute zero, since only con- 
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STATISTICAL THEORY AND THERMODYNAMIC 
FUNCTIONS 

A rigorous development of statistical mechanics requires a 
treatise of its own, and will not be attempted liere.f Elemen- 
tary statistical theory is, however, required in order to 
follow (1) the thermodynamical theory of the solid state, 
especially the heat theorem, and (2) the computation of 
thermodynamical functions for polyatomic gases from spectro- 
scopically determined energy levels. 

In applying statistical theory to the calculation of thermo- 
dynamic functions two main issues arise. In the first place, 
it is necessary to know all the energy states of any impor- 
tance which can be 'occupied 5 by the molecules in the system 
under discussion. This is a physical problem, and does not 
require statistics. As will be seen in what follows, for solid 
systems it is necessary to devise suitable approximations 
to these energy levels in order to calculate the thermo- 
dynamic functions at all. For gases, on the other hand, the 
object is to use the experimentally determined energy levels 
from band spectra. The second and purely statistical problem 
is to determine the average number of the molecules of the 
system occupying any energy state. Once this average 
number is known it is fairly easy to calculate all the thermo- 
dynamic functions required. 

Boltzmann’s theorem for the most probable distribution 
of molecules 

In statistical theory it is much easier to calculate the most 
probable distribution of molecules amongst the different 
energy levels than the average distribution. The difference 
arises from the fact that the molecules of the system need 
f Cf. Fowler, Statistical Mechanics , 2nd ed. Cambridge, 1936. 
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not exist in the most probable distribution all the time, but 
might fluctuate and occupy less probable distributions part 
of the time, thus altering the average distribution. The 
theory of fluctuations shows, however, that such departures 
from the most probable distribution are insignificant pro- 
vided the number of molecules over which the average is 
taken is sufficiently large, i.e. that it is quite safe to calculate 
thermodynamic functions from the most probable distribu- 
tion, provided the calculation is restricted to macroscopic 
systems. 

An elementary proof of Boltzmann’s theorem for the 
most probable distribution of molecules among different 
energy levels may be developed as follows: 

Consider a system with possible energy states 
€ 0 ’ ^ 2 ?***> 


and with N molecules in all. Any one distribution may be 
characterized by saying that it refers to n 0 molecules in the 
state e 0 , n r in the state e r , etc. 

If there are no special restrictions (cf. p. 94), the number 
of ways of realizing this distribution is simply 


njln 2 ln 3 l ...n r l ..." 


The most probable distribution is that for which this expres- 
sion is a maximum, i.e. the probability of a distribution is 
assumed proportional to the number of ways of realizing it. 

Any distribution is, however, subject to a number of 
physical conditions. In the present case these are 

(1) That the total number of molecules is constant, i.e. 

2 n r = N = const. (2) 

(2) That the total energy is constant, i.e. 

^ n r e r = E = const. (3) 

In order to determine the maximum probability W e of W nr 
it is necessary to have BW nr = 0 , or 

3 In W nr = 0 = 3 In A! — 2 Sin n r L 

H 


4360 


( 4 ) 
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When the numbers n r are large, the factorials may he 
evaluated, using Stirling’s theorem, in the form 

1 nn r l = n r lnn r — n r . (5) 

When they are not large, statistical theory either does not 
apply in this simple form, or the states in question make a 
negligible contribution to the behaviour of the whole. 

From equation (2) 2 & n r = (2') 

and from (3) 2 * r $ n r = 0. (3') 

Applying (2') and (5) to (4), since N is constant, this gives 
Slnl^ r = 0 = — 2 §n r \nn r . (4') 

Multiplying (2') and (3') by arbitrary multipliers A and ju. 
and adding to (4'), these equations may be summarized in 

the form £ 8n r (lnn r +A+^e r ) = 0, (6) 

which includes all the restrictions on the variations 8n r 
applying to the present problem. Since the values of 8n r in 
equation (6) are quite arbitrary apart from such restrictions, 
the only way the equation can be satisfied is for each of 
the equations \nn r +X+^ r = 0 

to hold separately, i.e. 

n r — Be-P*', (7) 

where B = e -A . 

This equation gives the most probable distribution of 
molecules in a very simple case, but Boltzmann’s theorem 
applies to a very large number of physico-chemical problems, 
unless these involve other restrictions than those considered. 
In order to use it in practice two points may be mentioned: 

(a) Evaluation of the constants A and fju in terms of physical 
' magnitudes. For the example chosen this can most readily 
be done by comparing the expression obtained with the 
concrete case of a perfect gas. Thus the total energy E of 
the N molecules is 

E = ^n r € r = B 2 e r e-^, 
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and the average energy EjN per molecule is 


B^€ r e-P'r 

~b2 


= e. 


i.e. 


e 


E 2 € r e ~ lxer $ 

N ”* "2 e“^ er ~ 


In 2 


For a perfect monatomic gas the energy states e r refer to 
the translational energy, i.e. 


€ r = +^I +^1 ) , 

where m is the mass of the molecule -and p l3 p 2 , p 3 are the 
momenta corresponding to the rth state, measured along 
any three axes at right angles. Since these can take up a 
continuous range of values, the summation may he replaced 
by an integration, i.e. 

N = 2 n r = B 2 e~^ r 

= JJJJJJ e^ 2m) ^*+ p ** piz) dp 1 dp 2 dp 3 dx 1 dx 2 dx d 

= Bv J JJ dp 1 dp 2 dp z , 


where v is the volume containing the N molecules. 

If E is the total energy of the gas, the value of p for any 
single molecule must He within the limits 
^J(2mE) > p > — *J(2mE). 

These limits when inserted in the integral give values prac- 
tically the same as when the integration is carried out from 
+oo to — oo, owing to the shape of the curve y = e~^ 2m)pS , 
Inserting the limits oo, — oo, which give a definite value to 
the integral. 


Similarly, 


whence 




For a perfect gas, 


N = 2 *,(^r) # . 

v Bv 3 ('n % {2m) s \ i m 

~ 2m 2 \ M 5 j ’ 

= 3^ 

2 E’ v \2irm) ' 


PF = |JS? = JV/ M = RT 
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anCt ^ = TcT' 

Thus the most probable distribution of molecules among 
the different energy states is given by the expression 

n r = Be~ €y l kT , 

where B is a constant determined by the physical properties 
of the system. The fact that 

fJ ' = kT 

depends on the role of \jl in the various thermodynamic 
functions for any system of molecules, when these are 
expressed in terms of the Boltzmann distribution (cf. pp. 53 
seq.), and is true for every physico-chemical system for which 
this distribution holds. 

Quantum weights 

In a number of physical problems a group of energy states 
€ r , e r+1 ,..., e r+q can arise in which the energy difference 
between e r and e r+q is very small compared with the interval 
to the next group of states. The group of q states is then 
conveniently taken together as a single state of quantum 
weight q> i.e. instead of writing 

n r = Be~* rlkT } 
n r+ 1 — Be~ €r + llkT , etc. 
the expression n r = Bq r e~ CrlkT 

is used for the whole group, which for purposes of calculation 
is equivalent to a single state of quantum weight q r . 

The partition function (Planck’s ‘Zustandsumme’) 

In any physical problem the distribution of molecules in 
the different energy states can be conveniently summarized, 
in virtue of Boltzmann’s theorem, by the partition function 
Z = q 0 e~ € ^ kT +q 1 e ~ € ^ kT +.,.+q r e~ € rl kT ... . 

This expression merely sums the number of molecules in 
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each energy state, but as a result of its structure it can be 
directly connected with the thermodynamic functions already 
described. 

Particular attention should be paid to the different evalua- 
tions of the energy levels used by different authors. One 
consistent method is to use the free atoms in the ground states 
as zero (cf. p. 57 ), in order to express the energy of the 
molecules. As defined here, the energy € 0 of the ground 
state of the molecule is merely equal to the spectroscopic- 
ally calculated zero-point energy. The only practical 
advantage this has over using the free atoms is that the 
dissociation energy is not always known as accurately as the 
energy differences e r — e 0 . 

Relation between Z and the number of molecules N in the 
system. 

According to Boltzmann’s theorem the number of mole- 
cules in the rth energy state may be written 

n r = Bq r e- CrIkT , 

where B is a constant of proportionality. The total number 
of molecules in the system is thus N = FJ q r e~ €rlkT } i.e. 

N = BZ. 


Relation between Z and the internal energy E. 

The average energy e per molecule is given by the expression 

-E_ B 2 e r q r e-^T 

N~ BJ,q r e-^ kT • 

The numerator is , ; the denominator is Z, 

8{l/hT} 8 y. 

1 8Z ain Z 

16 ' 6 ~ Z 8[i dfj, ’ 


E — Ne — -N 


a I nZ 

8/u. 
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since 


For a perfect gas , H = E+RT per mole. It is clear that 
this is equivalent to writing 

H = -N — 

djj, 

dlnfji _ 1 _ RT 

dfx f.i N 

Relation between Z and the specific heats. 

BE 0/x __ BE/ 2 BE 

dfji 8T~ 8y,\ TcT 2 ' h Bfx 

„ 0 B 2 lnZ , 

s= Ru? — 0 - per mole. 

OfX 1 

n BH ^ 2 B 2 l xi(Z/ijl) , 

Gp = 8T = ^ P er mole - 


r _8E 

Cv ~8T 


Complete partition function for several modes of energy intake. 

When the molecule can take up energy as translational, 
vibrational, rotational, and electronic energy, these different 
energy intakes can frequently be treated as independent to 
a first approximation. When this is so, the complete parti- 
tion function 

Q- = Z^ ran s X Z v jk X -^rot X ^elect X . • . 
is defined as the product of the separate partition functions. 
The reason for this definition is at once clear when it is 
noticed that the separate partition functions always enter 
into the expressions for the thermodynamic functions in the 
form In Z\ if the different energy intakes are independent, 
they must contribute in an additive way to functions such 
as E. This is secured by taking the product of the partition 
functions for G (cf. p. 116). 

Relation between Z and the available energy A. 

. From the Gibbs-Helmholtz equations (cf. p. 15), 

dA 
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i.e. 


E = A — T 


d_A t 

dT l 


' T 2+ TdT 


or, dividing by T 2 , 

_E_ _ 

T ' 2 ~ 

so that 

and using the value for 


E = — AT 


dlni5 

~dyT' 


we have 
or 


A = — i^lnZ. 


d_(A 
dT\T )' 

— k ^ E dfju; 


Proceeding in exactly the same way for j F, it is clear that 
for a perfect gas F = -RT ln(^//a) . 


From the property mentioned on p. 21 it follows that in a 
system in equilibrium Z must be the same for any component 
in every phase. 


Choice of zero in the evaluation of thermodynamic functions 

from Z. 

A problem of some delicacy is the choice of zero for the 
thermodynamic functions, when these have been evaluated 
from Z and are to be used in conjunction with thermo- 
dynamic functions evaluated in other ways, e.g. from thermal 
data and the heat theorem. Until the functions evaluated in 
different Ways have been adjusted so as to be measured from 
the same zero, it is convenient to refer to the thermodynamic 
functions evaluated from thermal data, and referring to the 
crystalline solid at absolute zero as standard, as H.T. func- 
tions (heat theorem), and to thermodynamic functions 
evaluated from the complete partition function as T.N.S. 
functions (theoretical functions including nuclear spin 
effects; cf. p. 104). 
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The difference between 
clear from a consideration 

HEAT 

CONTENT ENTROPY 


OBSERVED HEAT 
CONTENT - COOLING 

AS LOW AS POSSIBLE 
( EFFECTIVELY 0°K) 

APPARENT 
H.T. ZERO 

OBSERVED 
ENTR0PY<00LIN0 
AS LOW AS 
POSSIBLE 

TRUE 

H.T. ZERO 

UNDETECTED 

SPECIFIC 

HEAT 

ANOMALIES 

ZERO POINT 
ENERGY 

STATISTICAL 
ZERO (T.N.S.) 

ENTROPY DUE 

TO NUCLEAR 
SPIN & ISOTOPE 
MIXTURES 


Tig. 5 a. 


H.T. and T.N.S. functions will be 
of the different ways of evaluat- 
ing the internal energy and 
entropy of a gas, since the other 
functions can be made up from 
these two. The H.T. function 
E refers to all the energy which 
would be given up on cooling 
the gas down to absolute zero 
at constant volume. Inspection 
of Fig. 5 a shows that the T.N.S. 
function includes in addition 
the zero-point energy of the 
molecules, which is not thermo- 
dynamically observable in the 
mere process of cooling, i.e. 

^h.t. = ^t.n.s .-^° 5 


where E° is the zero-point energy per mole. 


Adjustment of zero for energy . 

In evaluating equilibrium constants from the equation 

—RTlnK = A = E-TS 

it is necessary to use the correct zero for both E and S . In 
thermodynamic processes involving only one molecular 
species the zero for E cancels. In a chemical reaction, the 
zero for each molecular species must be such that the alge- 
braic sum for products and reactants gives the thermo- 
dynamically observable heat of reaction at 0° K. With the 
partition function as defined on p. 53, this heat of reaction 
is automatically included in the equations on p. 58, but 
if A is evaluated calorimetrically or from Giauque’s parti- 
tion function (p. 60) it has to be introduced in the equa- 
tions on p. 58 by a zero-point factor e~ AJE ° !kT . (See note on 

p. 60.) 
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Adjustment of the zero for entropy in the complete partition 
functions. 

A farther complication, which must be adjusted before 
T.N.S. functions can be used in conjunction with H.T. func- 
tions, is that the complete partition function G includes 
factors which are never affected by thermodynamic pro- 
cesses such as change of temperature or concentration. One 
of these is the nuclear-spin multiplicity (p. 117), which is 
not changed in thermal or chemical processes affecting the 
atoms; in consequence, the factor cancels out in all thermo- 
dynamic calculations. Another is the isotope effect. Accord- 
ing to definition, the entropy of a mixture of two isotopes in 
proportions x : (1—x) would be 

S = ^^+(1 — x)S 2 — xhix — (1 — a;)ln(l — x), 

where x\nx-\-(l — rr)ln(l — x) is the mixture term, but since 
no ordinary process affects the ratio of the isotopes the 
mixture term usually cancels. 

If T.N.S. functions were used exclusively, no harm would 
be done by leaving in these factors in G , since they cancel 
whenever the functions are used in a calculation. They are, 
however, not included in the H.T. functions since e.g. the 
nuclear-spin multiplicity would not contribute to the specific 
heat till very low temperatures indeed. 

The safest procedure is to subtract all terms from the 
T.N.S. functions which refer to effects which could not be 
observed in ordinary ^thermodynamic processes (cf. also 
pp. 117, 120). 

It may be noted that some authors use the free atoms as 
zero of energy for the energy levels of the molecule. When 
this is done, all these levels differ from those discussed above 
by where is the dissociation energy of the molecule 
in its ground state, into unexcited atoms. The effect is 
simply to multiply the partition function by the term 

^atom zero ~ e ^molecule zero 

4360 T 
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and the adjustment of zero of the T.N.S. functions to corre- 
spond with H.T. functions offers no difficulty. In this case 
the zero-point energy factor is automatically included in the 
partition factor for vibrational energies. 


Expression for F in terms of G. 

After adjusting the zero of the complete partition function 
G by subtracting factors which refer to processes not observ- 
able thermodynamically, the complete expression for A F in a 
reaction a A +bB+... ^ mM + nN + ... is obtained as follows: 
For the individual substances, 


where 



includes a constant a, the partition function t for 
internal energies which are thermodynamically observable, 


and the factor 


(2TrmlcT\% 


for translational energy in three 


degrees of freedom. The constant a — 1 atm. in dynes/cm 2 . 
= 1*0132. 10 6 is inserted since F standard must refer to the gas 
at one atmosphere pressure. 

Also AF standard = — RT In 

so that 

p O a A G b B ..\a ) 


2 > 

-'Ll) 


By omitting the term 
the expression 


^/pj(Sn-S6) 


since K n — K. 


(Zn~Zb) 


Qm qu 

K c = na 7W~~ (concentration in moles/c.c.) 

is obtained. It may be noted that the unit of pressure in K p 
is the atmosphere, but that K p may be obtained in c.g.s. 
units by omitting the factor a. 
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Relation between the 'partition function and the entropy . 
Using the thermodynamic relation 


S = 


E—A 
T ’ 


it is easily seen that the value of S in terms of Z is 


a t>v r N din Z 


RlnZ — R 


din Z 
^ln/x 


A more elegant expression for S was first obtained by Boltz- 
mann, viz. 5 = k]nWej 


where S is the entropy of the system, and W e refers to the 
probability of the most probable distribution of molecules, 
which has already been discussed (cf. p. 49). This result 
may be obtained as follows : 

Using the expression on p. 49, 

lnT^ = Ain N— N— 2 ^hin r — J n r- 

But n r = J Be~ €rlkT , where the q r states of energy e r are all 
counted as separate, since they occur separately in W e . Also 
N = 2 n r = B 2 e^rlkT = JBZ. 

Substituting, 


]nW e = NlnZ+NlnB-B £ e-^ln Be ~ €rlkT ^ 
i.e. =N]nZ+2^. 


or k]nW e = = ^t.n.s. 


This expression makes the statistical significance of the 
heat theorem clear. If S -> 0 when T 0 this means that 
PTg -> 1 when T -> 0, i.e. that all the molecules of the crystal 
pass into a unique energy state at the absolute zero. 

At present it is not possible to determine experimentally 
whether the heat theorem applies only to or also to 

i* e * whether the minute energy differences leading to 
liuclear-spin multiplicity would be resolved at 0° K. 
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In passing, it may be noted that the various formulae just 
given make certain peculiarities of the thermodynamic func- 
tions much clearer. 

Thus the fundamental formula 
S = kh. W e 

corresponds with the fact that entropies are additive and 
probabilities multiplicative. When a system moves into a 
more probable state, its entropy is increased and its proba- 
bility is multiplied by a factor. 

In the same way the equation 

A = E-TS = -BTlaK 

corresponds with the fact that equilibrium depends both on 
the energy and probability of a system. At 0°K. the state 
of equilibrium is that for which A" is a minimum, since 
S = 0, but at higher temperatures states of bigger proba- 
bility and hence bigger entropy can be favoured in spite of 
the fact that they also have a bigger energy. 

A further variation in writing the partition function is used e.g. by 
Giauque (J.A.C.S. 52, 4808 (1930)), in which the zero-point energy is 
not included in the energy levels, so that 

z o = 1o+W* l/,cT +W- e,lkT +-- 

If it is desired to use this function the change can be readily made 

z,x«-jw = z . 

in all the formulae quoted. * 
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OwrsrG to their important role in the heat theorem, theo- 
retical expressions for the thermodynamic functions of a 
solid have received considerable attention. Apart from 
special effects, which will be mentioned later, the energy 
content of a solid is due to the vibrations of the atoms 
and molecules in the crystal, and the problem is to find 
workable approximations for the energy states which the 
atoms can occupy. Once such approximations have been 
found the distribution of the atoms amongst them follows 
the principles already discussed. 

A. Thermodynamic functions of a single oscillator . 

The simplest type of vibrational energy is that of a set 
of Planck oscillators of constant frequency v , and capable of 
taking up quanta of energy hv, 2hv,..., rhv and so on. 

For such a case 

Z = g 0 e~^ kT +q 1 e~ (<r ° + hv ^-{- „.+q r e - ieo+rhv)lkT + . . . . 

All the quantum weights are unity, and the resulting infinite 
geometrical progression may be summed to 
Z = e -hvlkT ) 3 

giving e = hvj(e hvlkT ~- 1 ) + ihv. 

The term hhv refers to the zero-point energy of the oscillators. 
The available energy per oscillator: 

a = kTlnil-e-W^+ihis. 

Applying these results to crystalline solids, if all the atoms 
in the solid vibrated independently with the same frequency 
v , the energy content of N atoms would be simply 
E == Ne = Nhvl(eWW-l)+\Nhv. 

This elementary theory was first proposed by Einstein, but 
although it gives a specific heat G v = BEjcT falling to zero 
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towards absolute zero, as required by the heat theorem, the 
rate of decrease of the Einstein function is more rapid than 
experimental results actually show. The failure of the theory 
is due to the fact that the vibrations of the atoms in a solid 
lattice are not independent but coupled, since the displace- 
ment of any atom from its position of equilibrium affects the 
restoring force on each of its neighbours. For some purposes 
the thermodynamic functions of an ‘Einstein’ solid may be 
used to discuss its behaviour, especially as the detailed 
calculation of the frequencies arising in a three-dimensional 
solid is very laborious. 

B. The ‘line’ solid with coupled vibrations. 

It is frequently useful to consider the simpler case of a 
one-dimensional solid. Taking a row of atoms of equal 
masses (the so-called line solid) and assuming that the restor- 
ing force on any atom when displaced from equilibrium is 
solely due to the displacement relative to its two immediate 
neighbours, for a simple harmonic vibration of frequency v 
the restoring force on the rth atom, if the displacement is 
X r , will be 

4t rVI r = k{X r _ x -X r )+k{X r+ x -X r ) 

according to the above assumptions, since X r „ x — X r and 
X r+x ~— X r are the displacements relative to the two immediate 
neighbours. 

Writing as solution X r — a cos rC and substituting 
4i r 2 mv 2 a cos rC 

= aft[cos(r — 1)U — cos rO-j- cos (r-f-l)C — cosr<7] 

= 2aJbcosrC[cos C— 1], 

SI i 27T 2 ‘V 2 m 

i.e. cos C — 1 = 2 sin 2 JO — — ~ — , 

sin iC-wJfg), 

which gives a relation between the natural frequencies v of 
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the line solid and the values of C . These are subject to the 
further condition that, since the ends of the solid must be 
free from constraint, 




cos(n+l)G = 0, 


(n >j> 0), 


i.e. {n+l)C = i(2i+l)7r, C = ^±L tt 

2(n+l) 

and an analogous solution, 

X r = asinrC, 

gives 

-X,i+i = Sin(«+ 1)0 = 0, C- 2nj 

where j is an integer. 

If b is the spacing of the atoms, the wave-length (i.e. the 
distances at which the disturbance is again in the same phase) 
is given by ^ 

X= — b. 


(n+1) 


(n>j> 0),f 


from which the velocity of the disturbance, v = vA, 

= (1) 

The line solid shows the characteristic phenomenon of dis- 
persion, i.e. alteration of velocity with wave-length. For 
small values of C, A is very large and (sin C)J C c=L 1, 

v - = VS) from (i) - 

The smallest value of A is 2b, C — tt, for which 



showing that the long waves have the larger velocity. 


Three-dimensional lattice. 

Although a similar dispersion must hold for the three- 
dimensional solid, the calculation of the individual fre- 
quencies for use in the partition function is not usually 
practicable. 

•j* C — 0 is excluded since it makes the equations on p. 62 indeterminate. 
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A very successful first approximation to the problem has 
been developed by Debye, who treats the solid as an elastic 
continuum, i.e. neglects dispersion effects. For such a con- 
tinuum the number of stationary oscillations or natural 
frequencies between A and A+^A is 

dZ x = ^4t7 dX, 


where V is the volume of the continuum. 

Now if dispersion is neglected, i.e. taking 
v for the vibrations, 


v 

A’ 


dv = — v 


dX 

A 2 ’ 


a single velocity 


dZ v = (a) 

where v t are the two translational velocities of sound (taken 
as equal in an isotropic medium) and is the velocity of any 
longitudinal vibration. Taking a mean velocity v m such that 


3 _ _2 1 


dZ v 


12rrVv 2 dv 


In order to satisfy the further physical condition that the 
total number of degrees of freedom of a solid composed of N 
atoms is 3N, an upper limit i/ max is imposed on the possible 
frequencies such that 


3 N = 


12ttV r 
v 3 J 

o 


v 2 dv, 


or 


3 N = 


4t tV 


V 


3 

m* 


This equation makes it possible to eliminate the velocity 
and to write 9Nv*dv 


Using this approximation for the possible energy states of 
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the atoms in the lattice, and the average energy e v of an oscil- 
lator of frequency v , at temperature T (cf. p. 61), viz. 

hv , i r 

' = —i+i hv > 


the total energy content of the solid due to the vibrations 
of the lattice is 

Vm 

[iUr = J e v dZ v 

o 

9iY f" z , j hv \ 9 N fhv 3 dv 

■= * J’ * + ^J— • 

0 0 

For purposes of calculation this expression is transformed 
by substituting ^ hv __ 

IcT X 


to a generalized function of (®/T), i.e. 


e/r 

m,-" !y (D’J S5+^®- 

0 

In this expression the constant 0 is a characteristic 
temperature of the substance whose value depends on the 
elastic properties. This follows both from the direct deriva- 
tion, from which. ^ /fsN \ 

0 Jet/ \brV)’ 

or from the theory of dimensions, according to which 

m 


0 


A- 


Jc 3 


where M is the atomic weight, X the compressibility, p the 
density, and the constant of proportionality A is approxi- 
mately the same for substances of the same lattice structure. 
These expressions are occasionally useful in estimating the 
thermodynamic functions of a solid from calculated values 
of 0 in the absence of experimental data on the specific heat, 
but should be used with caution. Tables of the values of this 

4360 tt 
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integral for various values of (0/T) are given by Landolt 
and Bornstein.f The values of C v = dEjdT are obtained by 
direct differentiation of the expression for E, which gives 
(writing y = &/T) 


Or 



V 


J 


x z dx 
e x ~l 


2/ 4 

e y — 1 


'-0 -» 

At temperatures such that ®/T is small, the value of G v 
tends to 3i? per gramme atom, the Dulong and Petit value 
for the atomic heat of solids. At low temperatures, where 
®/T is large and oscillations of high frequency make negli- 
gible contribution to the specific heat, the integral 

eiT oo 


f x z dx 

J e x —l 


approximates to 


j 


x a dx it 4 

e x — 1 ~ 15' 


Thus 


^ = 15 


0 

7T 4 


and 


n 12 T 3 4 


aT z 3 


where a is a constant depending only on the characteristic 
temperature of the solid. 

It will be noted that this expression for the specific heat 
satisfies the postulate of the heat theorem that 


v dH 
lim -~ 
T-+0&T 


0. 


Por solids with simple cubic lattices the theory has been 
verified over a range of temperatures by calculating the 0 
value at each temperature from the observed value of C v , 
making use of the tables already referred to. If the theory 
were to apply rigorously, 0 should be independent of 
temperature. Actually a drift in the values of 0 is nearly 
always found, J though in many cases the change in 0 does 
not exceed 5 per cent, over practically the whole range 
where 0 varies with temperature. 

f Tabellen , 5th ed. % Cf. Blackman, Proc . Roy. Soc. 148 A, 365 (1935). 
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At low temperatures the result C v = aT z is found to 
apply even to complicated non-isotropic solids such as 
CaF 2 , (NH 4 )C1, SiC, 

and both crystalline and glassy quartz and glycerine. This 
is probably connected with the fact that at low temperatures 
the wave-lengths of vibrations which still contribute to C v 
are large compared with the lattice spacing, so that the treat- 
ment of the solid as a continuum involves fewer inaccuracies. 

Thermodynamic functions and equations of state for 
the Debye solid 

In view of the approximate nature of Debye’s theory, 
graphical integration of the experimentally determined speci- 
fic heat should be used in calculating the thermodynamic 
functions whenever suitable data are available. The equation 
of state of the Debye solid is, however, occasionally useful, 
and may be developed as follows: 

Assuming the heat theorem to apply to the solid, and 
choosing the solid at absolute zero as standard state, the 
available energy A at any temperature is given by 

A = -Tj^dT, 

where the value of E has to be inserted from Debye’s ex- 
pression e/T 

E = 9JRT f 

J e x ~~l 

(omitting zero-point energy, cf. p. 56). This gives 

v 

A = -E/3+3 BT j AL 
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where y = ®jT as before. This is the available energy of the 
solid at temperature T but still at the constant volume 
it occupied at absolute zero. In order to calculate the avail- 
able energy at temperature T and at zero external pressure, 
volume V , the solid may be compressed at 0° K. to the volume 
V, and then heated to T at this constant volume. If W is" 
the work of compression, 


A = W+3BT 



V z dy 
ev—l 



where the value of © (which varies with the volume) now 
corresponds to the volume V. 

To obtain the equation of state the equation 


P 


is used. This gives 
P = 


■a 


w_ dItT ±de 

dV T dV 


TI& 

W f yZ dy 

i©/ J ev-V 

0 

d& 


or multiplying through by V and noting that — = 

Co V 

on transposing, 


pV+V 


dW 

dV 


d In® 
'dlnV 

din© 

“dhTP 


9 BT 

= 


v 

IT\ 3 t 

\©j J 


YE r 


V 3 dy 
e v — 1 


© dln@ 

Fdlnr 


where E v is the thermal energy due to the oscillations of the 
lattice. From this equation of state such properties as the 

dV 

coefficient of expansion are readily obtained. 

It may further be noted that the change in 0 with V is due 
to the fact that when the solid expands the increased separa- 
tion of the atoms leads to changes in the mutual forces they 
exert on one another, and consequently to changes in 
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v = Furthermore, the Debye solid at low temperatures 
m h 


obeys the heat theorem, since A = 

C v = 12/3T 3 = ~T 


PT\ 

8 2 A 

dT r 


Other aspects of the vibrational heat content of solids 

Attempts have been made to improve on the theory of the 
solid treated as a continuum by allowing for the increase in 
velocity of the waves with increasing wave-length. Inspec- 
tion of equation (a) on p. 64 shows that in the 3 N vibrations 
available there must be relatively more of high frequency, 
where the velocity is smaller, than is allowed for in the theory 
of Debye. This probably explains the success of the empirical 
equation proposed by Nernst and Lindemann,’j' 

where E(v) is the energy/gm. mole for oscillators all with 
frequency v, i.e. (cf. Blackman, l.c.) 


E{v) = 


Nhv 

ghvjkT J 


. Nhv 
2 


(cf. p. 61). 


For solids made up of complicated molecules, such as ice, 
the vibrations of the lattice as a whole are often practically 
independent of the internal vibrations of the molecule. In 
such a case the Debye theory may be applied to the vibra- 
tions of the lattice and the Einstein expression to the 
internal oscillations of the molecules: \Ep\ is simply the sum 
of Debye and Einstein functions, 

1 [E v ] = Diej+E^+Ei ® 8 ), etc. 

At temperatures approaching the melting-point of the 
solid the displacements of the atoms from their mean 
positions no longer obey Hooke’s law, i.e. the restoring 
force is no longer proportional to the displacement. The 

t Ze.it. f. Electrochemie , 17, 817 (1911). 
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resulting anharmonic lattice vibrations increase the specific 
heat by a linear term 

Gy = 3jR-f -ocl 1 , 

though in practice it is seldom possible to disentangle this 
effect from other factors leading to an increase in the specific 
heat above the Dulong and Petit value. 


Other sources of the heat content of solids 


Besides being able to take up vibrational energy, the atoms 
and molecules in a crystal lattice can increase their energy 
content in a number of other ways. Such changes in energy 
of the atoms are of considerable importance for the theory 
of lattice structure and have a direct bearing on apparent 
departures from the third law of thermodynamics. 

For purposes of calculation, the simplest of such effects 
is that first discussed by Schottky. It refers to a solid in 
which the atoms or molecules can occupy two different energy 
states, numbered 1 and 2, with energy difference e. The value 
of e for any atom is assumed independent of the states of the 
neighbouring atoms, i.e. of the number already in the higher 
state. Even though such cases are seldom realized in practice, 
the calculation illustrates a number of features shown by 
nearly all energy transitions even when these are much more 
complicated. 

When the transition energy c is independent of the number 
N 2 of atoms in the higher state, the ratio relative to the 
number in the lower state is given by the Boltzmann 


expression 


§ = ie~‘ lkT > 


where q is here the ratio of the quantum weights of the two 
states. The average energy per molecule due to this effect 


will be 


N 2 € _ eqe-< kT 

N ± -{~N 2 14 -qe~ e l kT ’ 
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and the additional specific heat C vs due to this cause will be 
_ dN 2 _ € 2 e^ kT 

° vs ~ € dT ~ q hT 2 (1 +qe~^ kT Y 
per molecule. 

A plot of this additional specific heat is shown in the 
diagram. As for all additional specific heats due to energy 



transitions in the solid, a maximum value of C vs is observed 
in the neighbourhood of the characteristic temperature 
0 ejk and C VJ3 falls off to zero either above or below this 
temperature. 

Direct integration or an a priori calculation shows that 
the maximum contribution of such an effect to the energy 
content E and heat content H is simply Nq€/(l+q) per 
gramme molecule, since when T ;> © 

= qA\. 

When such transitions take place at temperatures below 
the lowest at which measurements of the specific heat have 
been made, the value of the heat content as calculated from 
specific-heat measurements will be too low by 

Nqef(l+q). 
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This would mean, however, that the characteristic tempera- 
ture © = e/Jc is considerably below the lowest temperature 
of measurement, and consequently that e is very small. In 
consequence (cf. p. 13) no serious error can arise from this 
cause in the evaluation of heat contents. 

Quite a different situation arises for other thermodynamic 
functions such as the entropy. The additional entropy due 
to the Schottky effect may be calculated from the integral 

AS = J C vs d In T, 

and if the integration is performed from 0 to oo (in practice 
this simply means over the range of temperatures where 
C vs is appreciable, cf. p. 71), this integral has the finite value 
AS = J81n(g+I). 

However small the value of e separating the two states, 
and, consequently, however low the temperature at which the 
transition takes place, if the effect has been missed through 
not carrying specific-heat measurements to sufficiently low 
temperatures, the calculated value of the entropy from heat 
data will be too low by 

J?ln2, jRln3, ..., etc., 
according as q = 1,2,..., etc. 

This calculation of the simplest type of transition in a 
solid may be compared with experiment in two ways. The 
additional specific heat will be superimposed on the normal 
curve due to the lattice vibrations. For cubical lattices the 
normal curve fits a Debye function with constant © with 
sufficient accuracy to permit a calculation of C vs from the 
experimental data, with’ the following results: 


Substance 

€ cals . ( per mole ) 

© {Schottky) 

© {Debye) 

Diamond 

2,120 

1,070 

1,840 

Silicon 

490 

246 


Grey tin 

137 

69 


Lithium 

407 

205 


Sodium 

189 

95 

159 

Potassium 

117 

59 

99-5 
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An additional check on the theory is that since this transi- 
tion leads to no volume change (there being no effect on the 
neighbouring atoms when any one atom changes its energy 
state), the coefficient of expansion follows Gruneisen’s rule, 


dV 

dT 


= YK 0 C v , 


where G v is the Debye specific heat obtained after subtracting 
G vs from the experimental value ; K 0 is the compressibility. 

The soundness of using the theoretical Debye curve to 
compute an anomaly has been questioned, f but at tempera- 
tures where the Debye specific heat is small in any case, no 
serious error seems likely to arise. Thus crystalline gado- 
linium sulphate, (cf. p. 89) shows a Schottky effect; the 
specific heat of the solid obeys the T 3 law down to 7°K., 
but below this temperature the specific-heat curve begins 
to turn upwards till at 1*6° K. the specific heat is 500 times 
the value to be expected from the T 3 law. 

In this case a satisfactory physical explanation is available 
for the energy changes occurring. The Ga +++ ion is in an 8 S 7/2 
state, i.e. it has an eightfold degeneracy in the absence of 
superposed electric and magnetic fields. Interaction with the 
non-uniform electrostatic field of the crystal lattice removes 
the degeneracy. This means that the eight levels have slightly 
different energies, AE 0-52 cal. /mole. In consequence the 
higher levels are progressively emptied in the neighbourhood 
of the temperature given by RT = 0*52 cal. Solid (ortho- 
para) hydrogen (cf. p. 106) shows a similar effect in the 
specific heat; the effect corresponds with the presence of 
three different energy levels (AE c=L 7*5 cals.) in the solid. 

From the standpoint of thermodynamics, more important 
examples of such effects arise when the transition takes place 
below the lowest temperature of measurement, and lead to 
an apparent departure from the heat theorem. The following 


t Cf. Fowler, Bakerian Lecture 1935, Proc. Roy. Soc. 151 A, 1. 
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tablef gives the difference between the entropy for a number 
of gases calculated from vapour-pressure curves or chemical 
equilibria, and calculated from the integral 

J C P d In T, 

assuming the heat theorem applies to the solidified gas. 


Gas 

Difference 'per mole 

Notes 

n 2 

co 2 

CO 

n 2 o 

NO 

— 0-07^0-20 
0-32±0*27 
1*06±0*26 \ 
1*12±0*10 J 
0*90zb0*32 
0*75zh0*10 

Within experimental error. 
Dubious. 

Two different observers. 


To appreciate the significance of this table it should be 
noted that jKln2 = 1*4 approx., and that a positive differ- 
ence in the second column indicates that in computing 
entropy from heat data some contribution has been neglected 
(cf. p. 28), leading to an apparent departure from the heat 
theorem of the expected order of magnitude. The results 
suggest that the molecules CO, N 2 0, NO undergo some energy 
transition in the solid below the lowest temperature of 
measurement. In view of the fact that the atoms in any of 
these molecules are very similar, but not identical, it has 
been suggested that the molecules take up random positions 
such as CO, OC in the crystal at ordinary temperatures. At 
very low temperatures the energy difference between an 
arrangement such as 

C— O C— O C— O C— O 
C— O C— O C— O C— O 
O—O C— O C— O C— O 


and a random arrangement, e 

•g- 


CO 

OC 

CO 

CO 

OC 

OC 

OC 

CO 

CO 

CO 

CO 

OC 


t Clusius, Nature , 132 , 775 (1932). 
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is comparable with. IcT and leads to a transition to the lowest 
state, with corresponding entropy and energy decrease. 

The ammonium -chloride effect 

When the transition into a higher energy level is not 
independent of the neighbouring atoms, complications arise 



TEMPERATURE - 
Fig. 7. 


which make exact statistical calculation difficult. This is 
due to the fact that the energy difference between two states 
depends on the' number of molecules already in the higher 
state. Furthermore, owing to the altered interaction between 
neighbours, the transition is accompanied by a change in 
volume and a change in the infra-red absorption of the crystal. 

The diagram shows the additional specific heat super- 
imposed on the normal curve in such a case. Such effects are- 
shown by most ammonium salts in the neighbourhood of 
230° K. and by other substances such as FeO (185°K.), 
Fe 3 0 4 (115° K.), MnO (116°K.), Mn0 2 (93° K.), ^-butyric acid 
(221° K.), HBr (88°K.),HI (70° K. and 120°K.),CH 4 (20° K.), 
and also by hydrides such as SiH 4 ,HF,PH 3 ,H 2 S,butnot H 2 0.f 

f Landolt Bomstein, Tabetten, 5th ed. 
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One explanation (cf. Fowler, l.c.) which seems particularly 
applicable to hydrogen compounds is that below the transi- 
tion temperature the asymmetric ions are oriented in the 
crystal, but that when they acquire more kinetic energy of 
oscillation than the maximum energy of mutual orientation, 
they begin to rotate freely in the lattice. It is clear that the 
orienting force on an ion surrounded by neighbours already 
rotating is smaller than when these were themselves oriented, 
which explains why the energy of transition depends on the 
number already in the higher state. The specific-heat anomaly 
of such interactions is much steeper than for the Schottky 
effect. 

In support of the view that the transition involved is a 
change from orientation in the lattice to free rotation is 
the increased symmetry of the crystal above the transition 
point. The volume change is likewise due to a decreased 
interaction between the ions. Furthermore, the hypothesis 
is supported by the fact that the different ammonium salts 
show the effect at approximately the same temperature; 
other substances show transitions at temperatures which 
are lower the more symmetrical the molecule, and H 2 0, 
which is strongly coordinated, does not show the effect, 
whereas H 2 S does. At the same time any change of transi- 
tion energy, due to the effect of atoms already in the state 
of higher energy, will give similar specific-heat curves. Ex- 
amples are the change from ferromagnetism to paramagnet- 
ism, or from an ordered to a disordered alloy, discussed 
below. The physical change in the oxides of the transitional 
elements has not yet been elucidated, but is unlikely to be 
some form of free rotation replacing mutual orientation of 
the atoms. 

The effect of such specific-heat anomalies on the thermo- 
dynamic functions is similar to that of Schottky transitions, 
but not so amenable to calculation. 
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General theory of specific heats in the neighbourhood 
of transition points 

From the standpoint of formal thermodynamics a transi- 
tion point in a substance which is connected with a change 
of state may be represented as the intersection of the free- 
energy curves of the two states. A well-known transition is 



the case of melting. The diagram shows an intersection of 
the F , P, T surfaces for ice and water by the F, T plane 
corresponding with P = 1 atmosphere. 

The state of lower free energy is the stable form, and since 
the curves intersect at 0*0° C. there is a spontaneous transi- 
tion at this point from ice to water. 

Although the free energies per unit mass (or per mole if 
the same molar unit is used in the two phases) are the same 
at this equilibrium point, the heat contents and entropies 
in the two states are not the same, so that the transition 
from ice to water is accompanied by an absorption of heat, 
a so-called latent heat. 

Even in the case of melting this formal treatment of a 
transition point is not wholly adequate, as is shown by 
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experiments on the specific heat of ice just below the melting- 
point, f Fig. 9 shows the sharp rise observed, though the 
range of temperatures over which it takes place becomes 
narrower as the ice is progressively purified. This rise 
resembles the ammonium-chloride effect, but its statistical 



interpretation is not the same, since allowance must be made 
for the effect of impurities. 

If AT m is the total lowering of the melting-point due to an 
impurity of concentration G 0 , and if no solid solutions are 
formed, a lowering AT when fraction X is molten will corre- 
spond to a solution of concentration C 0 /X. If Raoult’s law is 
obeyed, AT is proportional to the concentration of impurity, 
and A T 1 

AT m X ■ 


If L is the latent heat 
partial melting is 


Q- 


of transformation, the heat used in 


LX = 


L 


ML 

AT’ 


f Dickinson and Osborne, Bull. Bur. Stand. 12, 69 (1915). 
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and the anomalous specific heat below the melting-point due 
to the presence of an impurity is 


AC V = 


dQ 

dT 


AT 

L - y ' L m T 

~ A y7 2 ’ 


since 


dAT 
dT “ * 


This simple expression agrees with the observation that 
the specific-heat curves become steeper as the ice is purified, 
and experiment leaves it undecided how sharp the transi- 
tion from a really pure solid to a liquid can be. A rise of 
specific heat before the actual melting-point would involve 
an abnormally rapid change of slope of the F, P , T surfaces 
in the neighbourhood of the transition point, but as long as 
there is a latent heat these surfaces must intersect at an 
angle. 

In addition to the phenomenon of melting a number of 
other transitions in the solid state are known which are only 
completed over a range of temperatures, such as the change 
ferromagnetism-paramagnetism or ordered-disordered alloy. 
Such changes could not be discussed in terms of free-energy 
surfaces without introducing an additional degree of freedom, 
such as the degree of spontaneous magnetization or of spon- 
taneous order of the solid, and a statistical treatment is 
required for proper interpretation of the phenomena. 


Disorder ed-ordered solid solutions. 

Reference has already been made to the fact (p. 47) that 
mixtures apparently do not obey the heat theorem, but that 
certain solid solutions can decrease their entropy by ordering 
the atoms in the lattice. This effect has been observed for 
such alloys as brass, Ee 3 Al, Cu 3 Pd, Cu 3 Pt, CuAu, etc. X-ray 
investigation shows that above a certain temperature both 
the pure metals and solid solutions have face-centred cubic 
crystal lattices. On cooling the alloys sufficiently slowly, a 
fresh 'superstructure 5 appears below a transition range of 
temperatures, indicating that a haphazard arrangement of 
atoms on points in the cubic lattice has been replaced by an 



80 STATISTICAL ACCOUNT OF 

ordered arrangement in which all the atoms of one kind 
occupy corresponding positions, e.g. all the Cu atoms at the 
corners and the Au atoms at the centres of the faces. Order- 
ing can also he followed from other changes in the physical 
properties. For example, the electrical resistance decreases 
when a superstructure is formed in solid solution. 

From the standpoint of thermodynamics and statistics 
the process of ordering leads to a gain in energy, and tends 
to take place spontaneously on cooling. It is, however, 
opposed by the thermal motions of the atoms since it in- 
volves a decrease in entropy, and can therefore only become 
prominent at temperatures where the gain in energy on 
ordering is large compared with the average thermal energy 
of the atoms. It is only for certain alloys that these tem- 
peratures are sufficiently high for the requisite migrations to 
take place in a measurable interval of time. 

A simple formal treatment has been given by Bragg and 
Williams. f It is first necessary to define the degree of order 
in a solid solution. This can be done in a number of ways, 
one of which involves the idea of ‘replacements’. Starting 
with a cubic lattice composed entirely of B atoms, the solid 
solution involves the replacement of a certain number n of 
these by A atoms. If a fraction r of the replacements occur 
in the positions required by a superstructure for A atoms, 
there will be nr positions of order oc and n{l—r) positions of 
disorder j3. 

The degree of order S is defined as follows: If p is the 
probability that an a position is occupied by a replacement 
of A for B , 

1 — r 

i.e. when p = 1, S — 1 (complete order), 

p = r, S — 0 (complete disorder). 

The increase in potential energy at constant temperature 
f Proc . Roy . Soc. 145A, 699 (1934). 
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when a replacement is moved from a ft to an a position will 
be V = F(£ 5 T), a function of the number of A atoms already 
in a positions. 

The distribution of atoms corresponds with nr (B -> A) 
replacements among n sites, there being nr positions a of 
order. At any moment the number of a positions available is 

nr — ' pnr , 

and the number of ft positions available is 
*tt(l — r) — (1 — jp)nr, 

so that applying Boltzmann’s expression 

£_ __ r (l -P) e vlkT . 

X — p (1 — 2r-\-pr) 

If r — e.g. for Fe 3 Al, for which only the cube centres are 
positions of order, and x == V jhT , 

8 = tanh 

In order to obtain a quantitative solution, V is assumed 
proportional to the degree of order, 

V = V 0 S 3 

to allow for the interaction of the neighbours, but is assumed 
independent of temperature for simplicity, since there is no 
great expansion of the lattice over the range of temperatures 
considered. 

The degree of order at equilibrium 8 e will be given by the 
intersection of the curves V — V 0 8, and 
S = tanh VjJcT, 

i.e. S e = tanh ~~ . 

This equation gives a rapid increase of S e in the neighbour- 
hood of a critical temperature 



above which 8 e 0. For Fe 3 Al 

T c = 833° K., V 0 = 0*29 volts. 

M 


4360 
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The energy changes accompanying changes in <S e may be 
estimated as follows: 

dp = (1— r) dS, 


dE = Vnr dp = F 0 8nr(l—r) dS, 
and just below the critical temperature an approximation 
SI = 3 (T c -T)/T c . 


for 8, is 


Differentiating, 


2 8 e dS e = 


3 dT 
T n 


Substituting for dS from the above equations, 

~ = AC V — Qnrk(l-r). 

The maximum, for r = of r(l — r) gives AC Vm&x = \nh. 
For Fe 3 Al there are in all n — \N replacements per N atoms 
of both metals, so that the anomalous specific heat may 
rise to as much as J of the Dulong and Petit value for the 
lattice vibrations. More detailed treatment shows that the 
dissipation of order above T c is less rapid than is indicated 
by the simple treatment just given, so that AC V falls off 
more slowly above T c than is indicated. 


Paramagnetism and ferromagnetism 

The statistical treatment of paramagnetism is discussed in 
the following chapter. Paramagnetic solids usually obey the 
Weiss-Curie law, according to which the susceptibility 
X - C/(T+ A), 

where C and A are constants. At temperatures where T 
becomes of the same order as A, X is no longer independent 
of the field strength §>, and Curie’s law ceases to hold. Both 
C v and the change in heat content on magnetization become 
abnormally large near the Curie point, below which the solid 
becomes ferromagnetic. Typical Curie points are 755° C. for 
iron, 588° C. for magnetite, 376° C. for nickel, though for 
weakly ferromagnetic substances the Curie point is not 
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observed till much lower temperatures. When cooling curves 
are plotted for paramagnetic substances, a change of slope 
is observed in the neighbourhood of the Curie point, owing 
to the change in C v , but there is no arrest corresponding with 
the evolution of latent heat at a fixed temperature. 

In the absence of a strictly quantitative theory of ferro- 
magnetism, a formal thermodynamic treatment may be 
given as follows. If $ is the magnetic field and I the 
intensity of magnetization, the work done on the substance 
in producing this magnetization is J <r> dl and the heat con- 
tent may be defined as 

H = E+PV—%1. 


If [E+PV] is independent of the state of magnetization, 
then taking §) and T as independent variables, 


/as\ _ (dl\ 

W P,T W ^ / Jp,-d 


whence the entropy change on magnetization is 


\8TJ 




d§>. 


For paramagnetic substances I — X§, and if the substance 
obeys Curie’s law, X = C IT, so that 



which is only appreciable at low temperatures, 

$> 2 C 


A S = — 


2T 2 ’ 


and the additional specific heat of paramagnetic substances is 



Ferromagnetics . 

Ferromagnetic substances below the Curie point show a 
spontaneous magnetization even for § = 0. 
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If I s is the saturation magnetization, the inner field, 
following Weiss, may he written 

H m = AI S , 

where A is a constant, and the energy of magnetization 

E m = \AI% 

per unit volume, or %(A/p)I% per gramme. 



Fig. 10. 


The additional energy of spontaneous magnetization leads 
to an additional specific heat 


ACU 


3T 


A r dl s 
p ls dT' 
die: 


Both experiment and theory show that ^ is largest near the 


Curie point, so that the maximum contribution to the 
specific heat is in this neighbourhood, and is in general bigger 
where I 3 is still large. Thus the maximum value of A C Vm is 
observed at the Curie point, but A C Vm falls off more slowly 
on the low-temperature side of the maximum. 

Fig. 10 shows the addition to the normal specific heat near 
the Curie point of nickel. 



VIII 

THE ATTAINMENT OF VERY LOW TEMPERATURES 

The energy changes of atoms in a solid are in general charac- 
terized by the interaction of neighbours. In phenomena such 
as ferromagnetism, ordered alloys, or transformations of the 
ammonium-chloride type this interaction makes it difficult 
to give a quantitative calculation of the energy changes 
involved. The paramagnetism of many solids only shows 
evidence of such interaction of neighbours at very low 
temperatures, and this makes magnetic methods of cooling 
of great importance. 


Unattainability of the absolute zero. 

Before discussing the methods which have been used for 
attaining very low temperatures it is useful to appreciate 
the significance of a theorem of Nernst’s, which states that 
the absolute zero is unattainable by any reversible process in 
finite dimensions. 

This can be demonstrated from a Carnot cycle. At a 
temperature AT above 0°K. the change A -> B is carried 
out reversibly and isotherm ally: 


A S A . 


B ___ 


F-H 

AT 


The system is then cooled adiabatically to 0°K., in which 
process AS = 0, and the reverse process B A carried out 
at absolute zero requires entropy change AS — 0 from the 
Nernst Heat Theorem. Finally the system is heated adiaba- 
tically to AT. 

If such a cycle were possible, a spontaneous entropy 
decrease AS a ^ b could be produced at AT without any other 
change, e.g. in the working substance. Since this is contrary 
to the second law, the conclusion is that if the heat theorem 
holds it must be impossible to cool a working substance 
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adiabatically to absolute zero, i.e. the absolute zero 
unattainable. 

Inspection of the equations 


AS = Gy 


AT 

T 


A Q 
T 


is 


shows that this theorem is in a sense trivial. Provided 
reversible processes are available which give appreciable 
entropy changes down to the lowest temperatures, it is 
always possible to produce quite large changes of A T/T, 
since C v is very small. Prom the thermodynamic aspect the 
change from 0*1 to 0*001° K. is as important as the change 
from 100° K. to 1°K., since it is the ratio AT/T that matters. 
It is always possible to attain temperatures which appear to 
be very small numerically, but the theorem of the unattain- 
ability of the absolute zero automatically comes into force 
when no further entropy changes can be found, so that the 
ratio AT jT which measures the cooling cannot be increased 
further. 


Magnetic processes of cooling. 

Every substance has a diamagnetic susceptibility inde- 
pendent of temperature. This can be shown in an elementary 
way by considering the electrons in closed orbits as equi- 
valent to electric currents, which set themselves so as to 
oppose the action of an applied magnetic field. The atoms 
of a limited number of substances have in addition permanent 
magnetic ' moments, due to uncompensated spin or orbital 
magnetism, and such substances have a paramagnetic sus- 
ceptibility dependent on temperature. 

Methods of cooling depending on the liquefaction of gases, 
i.e. on the entropy change on vaporization, begin to fail when 
the vapour pressure of the gas becomes so small that it 
becomes technically impossible to produce a sufficient 
entropy change in a finite time. The conditions under which 
the paramagnetic susceptibility of solids can then be used 
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for producing entropy changes down to very low tempera- 
tures will be clear from the following discussion. 

The change in entropy on magnetization is 

AS =f$?L dS (of ' p - 88) 


and the heat evolved is TAS. 

In the magnetic method of cooling the paramagnetic 
substance is magnetized, and the heat evolved is absorbed 
by the surroundings, using helium as a transport gas and 
liquid hydrogen to absorb the heat evolved. When tempera- 
ture equilibrium has been reached, the helium is pumped off, 
and on removing the magnetic field a corresponding entropy 
increase must take place in the working substance, causing 
absorption of heat and a fall in temperature. 

In order that the working substance may give large 
entropy changes down to the lowest available temperatures 


it must have^ large. 


Since I == X§, where X is the sus- 


ceptibility, and since most paramagnetic substances obey the 


Weiss- Curie law 

az; 

dT 


X(T+A) = C, 


$>C 

(T+ A) 2 ’ 


AS = 


$> 2 C 

2(T+ A) 2 * 


It is clear that the entropy change on magnetization will 
only increase with falling temperature so long as T ;> A, so 
that a suitable working substance must have A as small as 
possible. 


Statistical deduction of the Curie law , and significance of A. 

The statistical deduction of Curie’s law for the para- 
magnetic susceptibility of solids is based on the assumption 
that the solid contains elementary magnets which tend to 
orient themselves parallel to an applied magnetic field, but 
which are prevented from a complete parallelism by the 
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thermal motions of the atoms. If the elementary magnets 
have moments fx, their potential energy when they are 
inclined at an angle 9 to the applied field § is -~§^cos0 if 
there is no interaction between neighbours. According to 
Boltzmann’s theorem, the number of molecules falling in a 
solid angle dQ ^ d(ms d(f> 

with their axes inclined at 9 to the field will be propor- 
tional to e /Lt^COS0/fcT^(Q 5 


and the resulting susceptibility will be 
fi cos 9 per atom. 

The average magnetic susceptibility per atom will be 
- 2 eix ^ cosdIkT d€l fx cos 9 

^ 2 ef *«cos tf/fca’dQ 

It should be noted that this quantity is obtainable more 
compactly in terms of the partition function 

z = 2 ^ cos0/fcT ^, 

since the potential energy e per atom in a magnetic field is 
e = — and since e = — dlog Z/d(l/kT) = — §/Z. 

According to the quantum theory, only definite values of 
cos 9 can be assumed by the magnetic moments. Except at 
very low temperatures or very large values of S3, the average 
susceptibility per atom is the same if classical statistics is 
used, i.e. if the summation is replaced by an integration. If 
there are n atoms per unit volume, the susceptibility X per 
unit volume may therefore be written 


Writing 


X = 


e ix% 008 QlkTp cos Q gj n 0 dQd<f>^ 

J J e^ C0S ^/ fc27 sin 9 d6d<f> 



x 


IcT' 


This reduces to 


X = nix cothir + i , 

L x . 
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and for small values of x gives the Curie law : 

X = VfL 

3kT' 

Before comparison with experiment the small diamagnetism 
of the solid must be allowed for. 

The Weiss correction only becomes appreciable when 
T A, and is due to the electrical and magnetic interaction 
of the elementary magnets with their neighbours. Suitable 
working substances for which A is very small are gadolinium 
sulphate Gd 2 (S0 4 ) 3 8H 2 0, iron ammonium alum 
Fe 2 (S0 4 ) 3 (NH 4 ) 2 S0 4 24H 2 0, 

and similar salts of the rare earths and the transitional 
elements. The interactions between neighbours are small in 
such cases because the electron orbits responsible for para- 
magnetism are deeply 'buried 5 in the molecule, and are 
surrounded by a very symmetrical electrical field from the 
neighbouring ions, owing to the large number of water 
molecules. 

From the thermodynamic point of view, it is interesting 
to discuss the source of the entrgpy of magnetization. In 
gadolinium sulphate the magnetic quantum number of the 
gadolinium ion is j = |, and in a very weak field the orbits 
of the electrons can take up angles given by 

COS 9=j/j, ( j—l)/j , ( i— 2)/i , (j— 3)fj 

with respect to the field. Allowing for the difference between 
clockwise and counter- clockwise spins, this gives 4x2 — 8 
positions. The magnetic moment per gramme atom is gjfi, 

Nek 

where /?, the Bohr magneton, = -j— — 5,567 units, and g = 2 

is the ratio of the mechanical to the magnetic moment of the 
electron. 

The work required to set the orbits parallel to the field is 
thus gjP$(l—cos6). If the salt is in the earth’s field only, 
this work is very small, and comparable with RT only below 

4360 -w- 
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10 _5 °K. Neglecting the electrical orienting forces, thermal 
agitation would lead to all states being equally occupied 
down to this temperature. (Electrical forces raise it to 
about 1°K.) If the specific heat could he measured down to 
this temperature, an anomaly would be observed whose 
integral { C v d InT = ~-Rhi8. This is missed in the ordi- 
nary exper im ental evaluation of entropy (cf. p. 28) since 
the temperatures are not carried low enough. In a strong 
magnetic field, on the other hand, the energy gain on setting 
the molecules parallel to the field becomes greater than RT, 
so that all positions are no longer equally filled. As a corol- 
lary, the entropy defect is now less than R In 8, the difference 
appearing as entropy of magnetization. This has thus a 
maximum value which cannot be exceeded, a result which is 
connected with the fact that the susceptibility begins to fall 
off as the solid reaches magnetic saturation (cf. p. 88), sp 
that the integral 

‘ _ ldl\ 

^magnetization - ® 


has a finite value. 


§i 



IX 

THE GASEOUS STATE 


Whereas progress in the thermodynamics of the solid state 
depends on the correct use of the heat theorem, the thermo- 
dynamic functions of gases have only become amenable to 
calculation by a systematic use of the theory of the partition 
function, together with energy levels deduced from band 
spectra. This is due to the fact that the experimental evalua- 
tion of even heat contents of gases at high temperatures is 
not an easy problem. 

The translational energy of gas molecules . 

The theory of the partition of translational energy among 
gas molecules is chiefly of importance in providing a theo- 
retical expression for the entropy constant of a gas. At 
ordinary temperatures where the equipartition of energy 
holds among the molecules, the translational specific heat of 
a perfect gas is 0 P = |i? per mole, and the entropy is 

s==mn^% 

where 91 is the entropy constant. 

If gases obey the heat theorem, i.e. have S = 0 at 0°K., 
the translational specific heat must fall off at very low 
temperatures. This phenomenon, which is called the de- 
generacy of gases, is very difficult to observe experimentally 
owing to the simultaneous importance of the van der Waals 
correction at low temperatures and high densities, at which 
degeneracy would be observable. 

Note on the van der Waals correction . 

From the experimental standpoint, when a gas shows 
considerable departures from the gas laws the results can 
be compactly summarized in terms of the fugacity function 
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This is defined to be such that the change in free energy on 
compressing a gas is accurately given by the equation 
&F = BTlnfJU 

where / l5 / 2 are the fugacities at P x and P 2 respectively. 

The correction to the value of A F which would be obtained 
if the gas laws held is thus 

A P = BT In P X /P 2 -BT In f x /f 2 . 

In the limit at very low densities the gas laws are obeyed, 
and / 2 = P 2 , so that the correction A P at any finite pressure 
is given by 

A P = PPlnP — BTlnf = BTlnP/f . 

Writing as the next approximation to the equation of state 
of a gas Ti/v 



BTlnf If' = BTln PjP' — a(P — P'), 
and in the limit (/' == P' 0) 

PPln/= BT In P — ocP, i.e. A P — ocP, 

Statistical calculation of degeneracy . * 

Although direct experimental observations are unsatis- 
factory, a theoretical calculation of degeneracy is possible 
on the basis of statistical theory. 

For the case of the most probable partition of energy 
among particles, two auxiliary conditions must always 
hold, viz. 

Constant total number of particles whatever the distribution : 

Zn r = N. 

Constant total energy : 2 n r € r = E. 

In order to define the energy states more closely, the 
molecules can be considered as located in a six-dimensional 
space by their coordinates, characterizing their distances 
from the origin, x l9 x 2 , x 3 , measured along three axes at right 
angles, and the projections of their momenta, p ± , p 2 , Pz, along 
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these axes. Any molecule whose coordinates lie between x x 
and x x +Ax x > etc., lies in a cell of volume 

Ax x Ax 2 Ax z A p x A p 2 A p 3 . 

According to any form of quantum statistics, these cells 
must have a minimum size A 3 , since by the uncertainty 
principle ApAx ^ A. From this restriction it is possible 
to calculate the number of cells corresponding to energies 
between e and e+Se. The volume of such a six-dimensional 
phase space is 

II I I I I dx i dx * dx z d Ps, 

and is limited by the fact that no molecule can he outside 
the physically available gas volume v, so that 


v = J J J dx x dx 2 dx 3 , 


and no molecule can possess energy greater than e if it is to 
lie within the phase space up to e , i.e. 

2me > pl+pl+pt 


This is the equation of a sphere of total volume 


V = ^-{2me)\ V) 
o 


which is the volume available from e = 0 to e = e. The 
volume available between € and e-j-Se is 


dV = 27Tv(2m ) f €*§€ . 


Since the minimum volume of any cell is A 3 , the number 
of cells a r lying between e r and € r +3e r , amongst which the 
molecules with this energy can be distributed, is 


a r 


dV 

A 3 


= 27 TV 


(2 m)* 

A 3 


e*3e. 


In classical statistics identical molecules are considered 
as distinguishable objects, for the purpose of counting the 
number of possible distributions, irrespective of whether 
they are in the same cell or not. The number of ways of 
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distributing the N r molecules in cells, n x in the first, n 2 in the 
second, and so on, is jy j 

n x ln 2 l ... ’ 

as has already been stated. 

Einstein-Bose or Eermi-Dirac statistics for degenerate gases 
both take into account a further limitation on the number 
of possible distributions of n r molecules of energy e r between 
the a r cells. 

According to the Einstein-Bose statistics, identical mole- 
cules can only be distinguished from one another if they are 
in diff erent cells. The term 'degenerate’, as in other cases, 
refers to the fact that it is impossible to distinguish two 
molecules of the same structure in the same cell. 

The number of distributions of n r molecules in the a r cells 
lying between e r and e r +8e r in this case will be 

(■ n r -\-a r — 1)1 
** 110 ,- 1 ) 1 ' 

This formula can be established in a number of ways. One 
of these, due to Brillouin, is as follows: 

Instead of considering the Einstein-Bose gas, consider the 
number of permutations of n r identical objects among a r — 1 
boxes, to be calculated as follows: The n r objects and a r — 1 
boxes are thrown haphazard along a line, giving {n r ~{-a r —l)\ 
permutations, or allowing for the fact that the n r objects are 
indistinguishable per se, 

(n r +a r —l)\/n r \ 

permutations in all. The last box is now added at the end 
of the line, and the molecules to the left of every box, and 
between it and its neighbour, are piled into it, leaving an 
arrangement of molecules in boxes which has 
(n r +a r — l)l/n r l 

permutations. This arrangement is clearly identical with the 
Einstein-Bose case, except that the order of the boxes along 
the line has no significance for statistics. When account is 
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taken of this order along the line, there are (a r — 1) ! permuta- 
tions of the a r — 1 boxes among themselves. The number of 
permutations of statistical significance is therefore 

(n r -\-a r — 1) ! 
n r \(a- 1)!> 

which is the required formula. 

This gives the number of ways of distributing the n r 
molecules in the a r cells with energies between e r and e r +$e r . 
The total probability of an arrangement with n x molecules of 
energies between e ± and ^+3^, n 2 molecules with energies 
between e 2 and <r 2 -(-Se 2 , and so on, is the continued product 
for all values of r : 


' _ TT (( n r+ a r 

n>r HU, Ma r - 


The maximum of W or log W gives the most probable dis- 


tribution. 


Using Stirling’s theorem (p. 50) and neglecting unity 
compared with a r , 

= 0 

— 8 [2 {n r +a r )\n.{n r -\-a r ) — n r lnn r — a r lna r ]. 

As before, the two auxiliary conditions are (cf. p. 92) 

8 2 n r = 0, 8 2 n r € r = 

and the distribution function in the Einstein-Bose statistics is 


_ a r 

%r ~~ Ce^r—i’ 

where C and /jl are constants to be determined as before by 
the physical values of 

N=^^n r , E = ]5> r e r . 

Fermi-Dirac statistics. 

The distribution just quoted may be compared with the 
classic Maxwell-Boltzmann formula 

M 

n ’ ~~ 

and -with the distribution obtained in Fermi-Dirac statistics. 



96 


THE GASEOUS STATE 


By analogy with Pauli’s principle for individual atoms 
the Fermi-Dirac statistics imposes the restriction on the 
distribution of n r molecules among the a r cells of energy e 
that each of the a r cells can contain one molecule or none 
and not the very large number (up to n r ) as in the Einstein- 
Bose statistics. In view of this restriction, the number of 
distributions of n r molecules in a r cells is simply the number 
of combinations of a r things taken n r at a time, i.e. 



Applying the same methods as before, the most probable 
distribution with log W = maximum gives 

n r = u r /(<7e^+l)* 


The entrojpy of a monatomic gas. 

The above expressions at once give the complete values 
for the entropy and internal energy of a monatomic gas, since 

N == 2 n r > E = ^n r <= r , 

and 8 = felnT^a*. 

For the Einstein-Bose statistics, substituting for n r and 
replacing the summation by an integration, since the number 
of terms is very large. 


N = A 


co 

J 




(Oe^— 1 )’ 


where A 


„ / 2m\ * 

- *”(w) ■ 


S — XbiO+E/T—A t «*ln(l— e-(“/C) *, 

0 

oo 

e* de 


. E = A 


r d 
J (Cei*- 


■ 1 )' 
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It may be noted that E is not quite the same as — n ^ since 

dfjb 

the Boltzmann distribution does not apply. 

The Fermi-Dirac statistics gives the similar expressions 

co 

N = A J e*/(Cer*+l) de, 

0 

CO 

8 = J e*ln(l+e-*“/C') <Ze, 

0 

oo 

E = .4 J e*/(Cef“+l) de. 

0 

Gas degeneracy will only become appreciable when these 
expressions differ notably from the Maxwell-Boltzmann dis- 
tribution, i.e. when ±1 becomes of the same order of 
magnitude as 

At temperatures sufficiently high to neglect ±1, i.e. where 
there is no degeneracy, and the gas laws hold, it becomes 
possible to calculate the entropy constant of a perfect gas, 
without having to decide which form of quantum statistics 
it would obey at lower temperatures. Thus 


N = (A/C) j eie-^de, 


whence 




M*T* 


X 3-074 XiO- 4 , 


where n is the concentration in moles/litre, and M is the 
atomic weight of the particles. 

Furthermore, where there is no degeneracy, i.e. CeV-* ;> 1, 
L t 1 \ e -#** 

l \ 1+ c^) ~ o ’ 


In 


8 = R\nC+EjT+(A/C) J e*e-/** de 

o 

= BTnC+E/T+R = JS(|-f-ln C), 
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since for a monatomic perfect gas E = -%RT. Substituting 
for C and writing p = nRT for a perfect gas, 

8 = toff+l+to — 

from which the entropy constant is seen to be 

§o _ ■, ( 2 ? Tmhfih 5 

R m h* ^ 

Thermal data of a degenerate gas; the specific heat of metals . 

The degeneracy of a gas becomes appreciable at tempera- 
tures where M 

C= - X 3-074 X10~ 4 
n 

becomes 1. As has been stated, for all ordinary gases the 
van der Waals effects at the same temperature are so 
marked that a verification of the statistical theory becomes 
almost impossible. It is, however, possible to discuss a 
number of properties of metals in terms of the Eermi-Dixac 
statistics. The atomic weight of the electrons is only 
M = 5-43 X 10" 4 . In problems of gaseous ionization (cf. p. 37) 
the electron gas is so dilute that the perfect gas laws apply, 
but the conduction electrons in metals are present at a much 
greater concentration. If every atom in the metal provides 
one conduction electron, the concentration n is of the order 
of 10" 1 moles/c.c., and in view of the very small atomic 
weight C 1 even at ordinary temperatures, making the 
electron gas in metals completely degenerate. The energy E 
of such a gas is given by the formula 
E = E 0 +*T 2 , 

where E 0 , the zero-point energy per mole' of the gas, is 

E = —l ^" Nh2 ‘ 

0 40 \ 7tv I m 

which for an electron gas is as much as 15 times the classical 
value %RT at room temperatures. The conduction electrons 



THE GASEOUS STATE 


99 


thus have very considerable mobility, but their contribution 
to the specific heat is only G Vt = a T = oEjdT, where 

« = ilpl-nbhK 
\9nJ ¥ 

This would become comparable with the specific heat due 
to lattice vibrations at a sufficiently high temperature, but 
since these lattice vibrations become anharmonic at higher 
temperatures and also contribute a term A C Vmk6 = aT to 
the Dulong and Petit value, it is not always possible to 
obtain distinct evidence of the specific heat of the conduction 
electrons at high temperatures. At low temperatures, on the 
other hand, the value of C v due to lattice vibrations falls off 
at a rate proportional to T 3 , whereas the specific heat due 
to the conduction electrons only falls off proportional to T. 
It should therefore be possible to obtain evidence of the 
specific heat of the conduction electrons by investigating 
metals at very low temperatures. Although the treatment 
of the conduction electrons as a degenerate gas is not 
strictly accurate, since it neglects the effect of the positive 
ions of the metal lattice, the specific heats of metals such as 
Ag at very low temperatures show the expected T term 
(cf. p. 27). 



X 


THE THERMODYNAMICS OF HYDROGEN AND 
THE ROTATIONAL ENERGY OF GASES 


The distribution of rotational energy in a gas can be dis- 
cussed in terms of the theory of the partition function and 
the energy levels deduced from band spectra. The thermo- 
dynamic functions of hydrogen and deuterium are, however, 
so dependent on a correct evaluation of the rotational energy 
that it is preferable to discuss these gases in detail before 
dealing with the special issues that arise for more compli- 
cated molecules. 

According to the classical theory of equipartition of energy, 
each of the rotational degrees of freedom of a molecule makes 
a contribution i-B to the specific heat. For polyatomic 
molecules whose rotation can be described about three axes 
at right angles the rotational specific heat is thus O rot = fjft. 

According to quantum theory, rotational energy is quan- 
tized in such a way that the energy of the rth rotational 
quantum according to the Schrodinger equation is 

r(r+l)A 2 
87 t 2 I 




where I is the moment of inertia of the molecule about the 
axis about which rotation is being considered. The moment 
of inertia of relatively light molecules such as 

F 2 25-3 x 10- 40 c.g.s., 0 2 19*5 X 10” 40 c.g.s., 

N 2 13*8 xlO - 40 c.g.s. 

is, however, sufficiently large to make the quanta of rota- 
tional energy and the characteristic temperature 

A _ A a 
rot 8tt*I1c 

very small. Equipartition is therefore set up for the rota- 
tional energy of molecules, with two exceptions. 
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In the first place, diatomic molecules and linear poly- 
atomic molecules have a very small moment of inertia and 
very large rotational quantum about the line joining the 
atoms, so that this degree of freedom does not contribute 
to the rotational specific heat except at very high tempera- 
tures. Such molecules therefore have only two rotational 
degrees of freedom taking part in the equipartition, and 

c I0 t = R- 

For simple molecules the effect of rotational energy is to 


add the term - to the chemical constant (cf. p. 37), 

tit S 


provided the temperatures are sufficiently low for the 
moment of inertia of the molecule to be constant. At higher 
temperatures the vibration of polyatomic molecules alters 
their linear dimensions, and moment of inertia, thus intro- 
ducing a correction to the rotational energy which cannot 
be expressed in any simple form. The only satisfactory way 
to calculate the contribution of rotation to the thermo- 


dynamic functions when this happens is by using the indi- 
vidual energy levels as obtained from band spectra (cf. 
p. 115). 


The thermodynamics of hydrogen 

The moment of inertia of the hydrogen molecule is only 
0*47 x 10 -40 c.g.s. units, and this makes the quantum of 
rotational energy so large that departure from equipartition 
and a falling off from the rotational specific heat begins 
to be appreciable not much below room temperature, and 
below 50° EL hydrogen has zero rotational specific heat, 
thus behaving like a monatomic gas. 

Although the magnitude of the rotational quantum was 
known, no satisfactory partition function giving a correct 
specific heat could at first be constructed. A diatomic mole- 
cule with rotational energy e r can have the rotational quanta 
distributed in 2r+ 1 ways amongst its two axes (cf. p. 114), 
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and thus has a quantum weight (2r+l). The partition 
function would thus appear to be 


where 

Writing 


^rot = f (2r+l)e~^, 
0 

__ r(r-\-l)h 2 

€r 8^I~‘ 

^2 0 _ 

“ 8tt Uk ’ T 


= 2 {2r+l)e- r ^+ x > a , 

Nh 2 cZln.Z , 

— - — per mole. 

8 tt 2 I da * 


dE ra f i 

~ per mole 
dJL 


dE vot da 
da dT 


as is usual when the partition follows Boltzmann’s theorem. 

Erom the band spectrum of hydrogen, 

I — 0*47 X 10~ 40 c.g.s. units 

, 85-90 

and or = — jp-. 

Like other thermodynamic quantities (e.g. the Schottky 
effect, p. 71) this expression for C TOt gives a decrease for C rot 
at low temperatures, with a peak in the neighbourhood of 
the characteristic temperature 0 = 85-90° K. This is quite 
unlike the experimental curve. 

The discovery of ortho- and para-hydrogen has shown the 
discrepancy to be due to the assumption that all transitions 
between the rotational states are possible under ordinary 
experimental conditions. 

If the nuclei have a spin, i.e. a magnetic moment with a 
definite direction in space, it is possible in symmetrical 
diatomic molecules (cf. also p. 118) to distinguish between 
molecules with spins parallel (ortho) and antiparallel (para). 
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From Schro dinger’s equation and an extension of Pauli’s 
principle it follows that unexcited H 2 molecules can only exist 
in the following states: 


Antisymm. Proper function. 
Nucl. Spins Parallel 
Ortho 

Odd Rot. States 1, 3, 5 


Symm. Proper function. 
Nuel. Spins Antiparallel 
Para 

Even Rot. States 0, 2, 4 


The quantum weight of the ortho state is three times that 
of the para state. It was first suggested by Dennison that 
if the transition H 2 ortho -> para was infrequent during an 
ordinary specific-heat determination, the correct energy func- 
tion to use is not 

E mt = 2 (2r+l)e r e-^r, 

but E I0t = £ 2 (2r+l)e r e- e 'l kT +% 2 (2r+l)e r e~ e '/ kT , 

r even r odd 

where the equilibrium ratio at room temperatures and above, 
ortho/para =3:1, is frozen in at low temperatures owing to 
the infrequent transitions, i.e. 

^rot ~ i ( ara ^ ^ortho ) P er mole, 

where ™ 

Cpaxa = Bo*^]n(I+5e-«<'+9e-™°+...), 


✓72 

^ortho = ^£- 2 ln(3e-^+7e-^ + ...). 

This equation gives excellent agreement with experiment, 
and can be further checked by working with different 
mixtures of hydrogen, containing from 100 to 25 per cent, 
para-hydrogen. In each case the appropriate <7 rot curve can 
be followed experimentally. It is interesting to note that 
for para-rich mixtures C TOt passes through a maximum 

^rot ^ 

Before discussing the effect of this 'frozen’ equilibrium on 
the experimental evaluation of the entropy constant of 
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hydrogen, it is interesting to discuss what causes a transition 
from the ortho to the para state. At low temperatures the 
change takes place slowly at high pressures, and with a 
catalyst of active charcoal the change to the appropriate 
mixture is complete within half an hour at temperatures of 
liquid air and liquid hydrogen. By pumping off the equi- 
librium mixture at the temperature of liquid hydrogen, 
practically pure para-hydrogen can he obtained, and its 
properties studied. 

The transition ortho ^ para is also catalysed by para- 
magnetic substances such as oxygen, NO, N0 2 , and ions of 
transitional metals and the rare earth salts. At high tempera- 
tures catalysts such as platinum black have an activity in 
favouring the transition which is parallel with their activity 
in hydrogenation. Hydrogen atoms also produce the change 

H 2 para + H ^ H 2 ortho + H 
Spins ^ -J4 t 

Differences in the physical properties of pure para-hydro- 
gen and the ordinary 25 per cent, mixture are small but 
definite. Thus the triple point pressure is 53*0±0*I mm. for 
para-hydrogen and 53*9±0'1 mm. for the normal mixture. 

Entropy constant of hydrogen 

The entropy of hydrogen presents certain difficulties of 
evaluation either when a purely theoretical method or a 
purely experimental method is used. The reason is that 
the theoretical calculation introduces terms for the contribu- 
tion of nuclear spin to the entropy, which arise in the same 
form for any change of phase of hydrogen, or for hydrides, 
but which are not included in any experimental evaluation 
of AS — J C r d In jP, since the energy transitions involved 
would not contribute to C P till a very low temperature 
indeed (cf. pp. 59, 72). 

In consequence, the theoretical entropy of hydrogen in- 
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eluding nuclear spin effects (T.N.S. entropy) gives correct 
values of equilibrium constants and vapour pressures when 
it is used in conjunction with the T.N.S. entropies of the 
other substances concerned, but not necessarily when it is 
used in conjunction with entropies calculated on the basis 
of the heat theorem and the specific heat integral (H.T. 
entropy). What happens is that the N.S. contribution 
cancels when T.N.S. entropies are used consistently, but 
cannot do so when H.T. entropies are used with T.N.S. 
entropies in one equation. 

A similar problem arises in a modified form whenever the 
thermodynamic functions calculated on the basis of thermo- 
dynamic laws (the H.T. functions) are combined with func- 
tions calculated solely from the complete partition function 
for all forms of energy, these being evaluated from spectro- 
scopic data (the T.N.S. functions). As a general rule the 
T.N.S. functions include effects which are missed in the H.T. 
functions, simply because there is no ordinary thermodynamic 
means of detecting them, i.e. simply because they do not 
affect the entropy changes in processes such as vaporization 
or chemical reaction, or changes in heat content (specific-heat 
determinations). In any particular problem it is always 
necessary to consider what portions of the T.N.S. function 
have to be discarded (this merely involves a change of zero of 
the T.N.S. function) before it can be used in conjunction 
with H.T. functions. Some examples are discussed below. 

The H.T. entropy of hydrogen. 

Since rotation makes no further contribution to the specific 
heat of hydrogen at low temperatures it can be treated as a 
monatomic gas. At rather higher temperatures the contribu- 
tion of the rotational specific heat to the entropy would be 
N rot = f C xot d\nT. For example, if the entropy of ‘mon- 
atomic’ para-hydrogen is taken as that of a monatomic gas 
at low temperatures (cf. p. 31), this would give values of the 
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entropy of vaporization of solid para-hydrogen (the only 
experimental check at low temperatures) in agreement with 
experiment. Furthermore, the rotational specific heat of pure 
para-hydrogen would give the correct H.T. rotational entropy 
at somewhat higher temperatures for any process in which the 
system contains only para-hydrogen. 

The H.T. entropy of the 3:1 ortho -para mixture at low 
temperatures does not refer to an equilihrium mixture, and 
must therefore he used with circumspection. If a process 
such as vaporization does not involve any change in the 
mixture , the experimentally determined entropy change 
would agree with the value determined from 

$3:1 ~ $monatomic“H ^rot d in ^ ' 

J 3:1 

Actually (cf. p. 73) the 3:1 mixture shows a Schottky effect 
in the solid state which has to be allowed for in calculating 
the entropy change on vaporization. Thus, neglecting the 
Schottky effect, the experimental entropy is only 29*65 units 
at 298° K., in lieu of 31*23. 

When the 3:1 mixture takes part in a process at very low 
temperatures, in which it is changed into the para form, it 
is doubtful whether thermodynamic functions applying to 
equilibrium phenomena can be used with advantage. The 
H.T. entropy is sometimes represented by 

$3:1 = l$ortho+i$para+fl n l + il n i» 

the last two are the mixture terms (cf. p. 46), but this 
expression has practically no connexion with experimental 
data. 

At temperatures between 50° K. and about 300° EL. where 
equipartition has not been completely set up for rotational 
energy, the correct H.T. function to use for the entropy is 

$ = i$para+l$ortho + iln| + llni, 

provided there is no change in the 3:1 ratio in the process 
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whose entropy change is being calculated. This also leads 
to a cancelling of the mixture terms. 

At temperatures above 300° K., C TOt = R } and hydrogen 
can be treated as a diatomic gas (at 298° K. the standard 
entropy is 31*23 units/mole). Applying the standard formula 
for the entropy constant (p. 37), 

8 = 22[ln2l+ln|+ln?^], 

s = 2, <7 = 27 + 1=1 (p. 115). 

The nuclear spin entropies may be added — using p. 1 1 8 — but 
are omitted here, apart from the symmetry factor (p. 120 ). 
The entropy constant of the mixture, i.e. of ordinary hydro- 
gen, is thus 

S = J^para + l^ortko + i In 1 In | 

= -ft £ln 9 I+ln — — — J-f-ln J+fhif+Jln 

where the last two "mixture 5 terms must be subtracted to give 
the H.T. entropy (cf. p. 57). ^ 

It is useful to note that the chief difficulty in evaluating 
the entropy of hydrogen arises at temperatures where the 
3:1 mixture is not strictly the equilibrium mixture, and that 
the calculated values can only be used for processes where 
the ratio ortho/para does not change. 

Thermodynamics of deuterium 

Two effects have to be discussed in connexion with the 
thermodynamics of deuterium: 

(а) A comparison of the rotational heats of H 2 , HD, 
and D 2 . 

( б ) A comparison of the chemical constants of H 2 , HD, 
and D 2 and their effect on the position of equilibrium. 
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The rotational specific heats . 

The problem for each gas is to construct the rotational 


partition function r? ^ „ 

"rot = z Qr e ~ 


■e r !kT 


since 


CL 


Rcr * a*inz TOt 

do 2 


The nuclei of hydrogen have spin i = J, of deuterium spin 
i = 1 (in units hj 2 tt) } so that both these gases exist in ortho 
and para modifications. HD is not a symmetrical molecule, 
so that rotational transitions from even to odd states are 
possible. The results can be summarized in a table: 


Gas 

h 2 

HD 


Nuclear spin. 

♦ = i 


i = 1 

Rotational states. 

para r — 0, 2, . 

all values 

ortho r = 0, 2, 4, 

Quantum wts. of 

ortho r = 1, 3, . 

r = 0,1, 2, ... 

para r = 1,3, 5, 

sy mm ./ant isymm . 
eigenfunctions 

1:3 


6:3 

Moments of inertia 
assuming con- 

stant separation 
of atoms. 

0-466 x 10" 40 

0*62 x 10" 40 

0-93 X 10~ 40 

Minimum quan- 
tum change. 

<N 

t 

O 

II 

r = 0 — > 1 

t 

O 

li 

£■» 


It will be noted that for HD the simple partition formula 
holds (cf. p. 102). This gives a sharp rise in the rotational 
specific heat in the neighbourhood of © rot and a subsequent 
drop after passing through a maximum, and this behaviour 
has been experimentally verified. 

Dor D 2 the rotational specific heat of the frozen equilibrium 
mixture (ortho/para 2:1) is 

^rot ~ §^ortho“^ s^para 5 

where C OTtho = Eo 2 -—In 2 (2r+l)e~ e ^, 

ao x even 

d 2 

^para •Kcr 2 -— In 2 (2r+l)e- e 'l kT , 

aa rodd 
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and 


for D«> 


^h 3 


h 


€ r for H 2 J Da 
The resulting curve is quite different from that for H 2 and 
has been verified experimentally, f 


Chemical constants of the gases H 2 , HD, and D 2 . 

The equilibrium constant of the reaction 
H 2 0+HD ^ HDO+H 2 

is 3-8 at 25° C. and 2*0 at 100° C. In general the equilibrium 
constant in an isotope exchange 

AH+BD ^ AD+BH, 

_ [AP][BH] 

[AH][BD] 

differs from unity on account of the following factors: 

(1) Differences in the zero-point energies of the isotopic 
compounds. Since ItTlmK — A F = AH — TAS, differences 
in the zero-point energies of e.g. AD and AH will affect 
In K by changing the value of AH when hydrogen is sub- 
stituted for deuterium. 

(2) Owing to the different masses of H and D, the entropy 
constant will differ. The magnitude of the effect will be clear 
from inspection of p. 36, and is most compactly formulated 
as follows: 

#d 3 = l^ortho+i^para+f 2 +^ 3 ( the mixture terms). 
The nuclear spin entropies are omitted apart from the sym- 
metry factor (p. 120), as are the mixture terms. 

• ** $d 2 = j^ln A+fln.-M'oa+h 1 r D2 j+ln|-, 

where A is a constant, 

= # H2 -f-j?[§ln 2]. 

Similarly, = ^H 2 +^[|ln|— InJ]. 

From the experimental standpoint, the effect of isotopes 
j* Clusius and Bartholome, Gottinger N achrichten, 1 (1934), 1. 
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on equilibria is readily studied by using hydrogen and 
deuterium and an exchange catalyst such as platinum. 
Equilibria in solution may be studied using H 2 0 and D 2 0 and 
acids as exchange catalysts. 

A reaction of fundamental importance is 


H 2 +D 2 ^ 2HD, 


_ [HD] 2 
2 [HJCDJ 


In the region of temperatures where the rotational motion 
is fully excited, and the vibrational motion is not yet excited, 


BT In K 


- ( 2e HD — e H, — «D») + f l 11 


Mr, 


+ ln J-|p . 

■Mb. 


2 

HD 


-In 4. 


The first term on the right-hand side represents the change 
in zero-point energies, which are 


€ = Tic(\v— Ivx), 

where v is the wave number of the vibrational frequency and 
X is the anharmonicity constant. The term In 4 arises since 
two symmetrical molecules are changed into two unsym- 
metrical ones. The experimental result is that 

2€ H d— ^h 2 — e D 2 = = — 155 cal. /mole 

and 1 nK x = -34/T+0-6276. 


The vibrational energy of gases 

One of the most general means of determining the specific 
heat of gases is from the velocity of sound in the gas. After 
allowing for corrections for departures from the gas laws, 
usually small under the conditions of experiment, 

C P -C V = E, 

and the ratio C P jC v = y is determined from the velocity 


u — 


where P is the pressure and p is the density. 
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After subtracting the translational and rotational specific 
heats, the vibrational specific heat of the gas can be estimated. 
For simple diatomic molecules only one mode of vibration 
is possible, and the partition function for vibrational energy 
might be expected to be the simple Einstein function (cf. 

P* 61) z = 2 e- 7iVo( - r+mT , 

r 

where \hv Q is the zero -point energy, and the quantum weight 
is assumed to be unity in the absence of complications. 

This result agrees with experiment only for the simplest 
gases, owing to a number of complications. One of these is 
that at high frequencies the velocity of sound has been found 
to vary with the frequency. For carbon dioxide at ordinary 
temperatures, for example, the velocity begins to increase 
in the neighbourhood of 10 5 cycles per second and finally 
reaches a new constant value. This is due to the fact that 
the vibrational specific heat takes a finite time to come to 
equilibrium distribution, which may become comparable 
with the time of an experiment in explosion methods, or for 
high-frequency sound waves. The vibrations of the molecule 
then no longer follow temperature changes during the 
adiabatic pressure changes. Since the velocity of sound is 
JR> 

aVy and y = 1+ — , the decrease in C v leads to an increase 
G v 

in the velocity of sound. 

It should be noted that a modified specific heat might 
also be more suitable in calculating phenomena such as the 
velocity of an explosion wave. 

A second reason why the simple partition function is un- 
satisfactory for the vibrational energy is that the individual 
vibrational quanta are fairly large (compared with rotational 
quanta, for example), and the vibrations become anharmonic 
when only a few quanta of energy have been taken up. This 
means that the value of the nth vibrational level is smaller 
than nhv 0 , and the values of the separate levels have to be 
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inserted in the partition function summation. These vibra- 
tional levels are deduced from the spectroscopic energy data 
for the molecule. 

The complete partition functions for some simple atoms and 
molecules . 

Apart from the s im ple theory already discussed, the diffi- 
culties in calculating the complete partition functions of even 
simple molecules are not due to thermodynamics at all but to 
the complicated nature of band spectra. 

The hydrogen atom. 

Eor the nucleus of an atom the partition function refers 
to the number of orientations of the nuclear spin in a per- 
turbing field. At all attainable temperatures the energy of 
orientation is so small compared with hT that each of the 
g~e r ikT terms in the usual partition function is unity, and 
Z 0l = 2i s +l, where i s is the number of units of nuclear 
spin (^/2tt). 

The only other energy intake to be considered is the elec- 
tronic energy. The quantum weight of any electronic energy 

lGVeliS It elect = (2j s +D, 

where j s = (0 3> j s )> 

and the electronic partition function is 

Select = 2 (2j s +l)e-*^. 

This function should be summed over all the electronic levels 
indicated by spectroscopic data, but actually e~^ s ^ kT makes 
an appreciable contribution to Z only when T > ej s l4Jc, so 
that only the electronic ground state is of importance, except 
where the separation of electronic energy levels is small. 

Eor hydrogen the electronic ground state has quantum 
number 1 and j 8 = The complete partition function (T.N.S.) 
for the atom is thus 

Z = Z nuc i X Z^ wf. “ 2x2 = 4. 
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For chlorine i s = § and i5 nncl = 6. 

The electronic ground state is an inverted doublet 2 jP. 
j s = f for the lower component. For the higher component 
Av = 881 cm. -1 and j s = For these two components of the 

doublet 1 = 4 and 2 respectively, 

so that ^eiect = 4+2e- 8817 ^ r 

and the complete partition function is 

% = ^nnclX^ el ect = 6(4+ 2e~ 88 *^). 

Quantum weights and degeneracy . 

There are several ways of illustrating the origin of the 
'multiplicity 5 of the principal energy levels of atoms and 
molecules. 

From the standpoint of formal wave mechanics, the multi- 
plicity q of a state with energy E arises when there are q 
distinct proper values satisfying the wave equation of the 
molecule in the particular quantum state. The proper func- 
tion of the molecule in this state is then a linear combination 
of functions with q arbitrary constants. The physical signi- 
ficance of degeneracy # may, however, be illustrated in other 
ways. 

The electronic quantum weights for atoms. 

The quantum weight of an electronic energy level is 
2J S + 1, where j Si the 'inner 5 quantum number, is the resultant 
of the magnetic moments due to the spins of the electrons 
and of the orbital magnetic moments. In vector notation 

j = s+1. 

The reason for inserting a quantum weight in statistical 
theory may be seen by considering what happens in a weak 
magnetic field. According to the principles of space quantiza- 
tion of magnetic energy, the resultant magnetic axis of the 
atom will set itself at such angles that its projection along 
the field can have the series of values 

m = j, j—l, —1, —f 

43<50 q 
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If the magnetic moment is the energy intake due to these 
orientations in a field of strength $ will be obtained from the 
usual partition function Z = ^ e^ m ^ kT , in which there are 

m 

Zj+l terms. In the absence of the magnetic field § = 0 and 
Z = 2j+l, which must still be inserted in statistical theory. 
Multiplicities may thus be regarded as due to the partition 
of atoms or molecules amongst the minor variations of a 
given energy level, and must be inserted in order to obtain 
the correct probabilities in the distribution function. 

The rotational degeneracy of diatomic molecules. 

A rotational level with energy E i = has the 

&7T 2 2 



multiplicity (2j+l), due to the (2?+l) arbitrary constants 
in the proper function describing that state. 

Another way of illustrating this arises from a considera- 
tion of the number of ways in which the j quanta can be 
partitioned amongst the two rotational axes of the diatomic 
molecule. The particular arrangements which can be chosen 
are shown in the diagram. When rotation is taking place 
about two axes at right angles, a distinction can be made 
between the case in which the rotation about the second 
axis is in the same or in the opposite sense. Its physical 
foundation is clear when the molecule is a dipole, since in 
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such a case the magnetic axes due to rotation can be dis- 
tinguished. 

There are^+ 1 distributions of they quanta between the two 
axes. Each of these, except that for y 2 — 0 5 must be doubled 
to allow for the sense of the rotations. The total number 
of minor variations of a rotational level is thus 

The complete partition function for some diatomic molecules . 

The summation has to be made over every possible 
electronic, vibrational, or rotational state contributing appre- 
ciably to Z. Since the ground electronic state is usually a 
1 S term, the multiplicity = 1. Common exceptions are NO, 
O^OH, CSKT. 

It may be noted that cases such as OH and CH are 
particularly interesting, since reliable estimates of their 
thermodynamic functions could hardly be obtained except 
from spectroscopic data. The thermodynamic functions of 
radicals are of particular value in verifying certain reaction 
mechanisms in which they are supposed to take part. Com- 
plications in applying the theory of the partition functions 
arise chiefly in the selection of the correct quantum weight 
for the various states. 

Approximate empirical formulae for band spectroscopic fre- 
quencies . ^ 

When the band spectrum of a molecule has been analysed, 
the results can usually be expressed in terms of an empirical 
formula 

v = v 0 +^ e ( v +i)~~ Xco e( v +i) 2 + — + 

+^J(/+l)+D r J 2 (J+l) 2 HWl) 3 +.., 

where 

vq is the electronic frequency separation, 
co e is the fundamental vibration frequency in any 
electronic state, 

x is the anharmonic vibration constant. 

The half-quantum of zero-point energy is added to v. v and J 
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are the vibrational and rotational quantum numbers. B v 
and D v vary with v, i.e. 

B v — B e — a(^+i)+'y( v +i") 2 *--5 
D v = A+^+i) 2 +S(^+i) 4 -... 

These formulae allow for changes of moment of inertia in 
the different rotational levels (the molecule may stretch under 
the centrifugal force due to rotation) and for the interaction 
of vibrational and rotational quanta (the average size of the 
molecule increases with v). 

Tor rigid molecules 

hr 

D v = E v = 0, B v = g^2/* 


Summation to obtain Z. 

Z is determined at a given temperature by summing for 
all rotational levels, keeping v = 0, and repeating for 
v = 1, 2, 3,... so long as the contribution to Z is still appre- 
ciable at the temperature in question. The process is re- 
peated for any electronic levels sufficiently close to the 
gro un d state to make an appreciable contribution to Z. 

In virtue of the exponential structure of the individual 
terms in the partition function, when the internal energies 
(electronic, vibrational, and rotational) can be treated as 
independent, the complete partition function factorizes into 
the factors 

Z = 2 2vib e~~ € * h!kT 2 ^rot e ~ €lot,kT 2 Select e ” %ect/fcr 
^ ^vib ^rot ^ elect? 


as has already been mentioned (p. 54). 

Integration methods can sometimes be used in place of the 
summation. Thus the simple rotational partition function 

oo 2 

2 (2J + l)e~ <,j ( jr+1 ), 


Z 


rot 


. 7=0 


8tt 2 I1cT’ 


can be evaluated by integration at relatively high tempera- 
tures where the interval between successive terms is small, 


^rot I/". 


givmg 
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It will be noted that this accounts for the rotational entropy 
term in the entropy constant for simple diatomic molecules 
(p. 37). 

The summation for the complete partition function must 
then be repeated over a range of temperatures, since the 
value of Z as a function of temperature is required for a 
complete evaluation of thermodynamic functions. 

Quantum weights in special cases. • 

In view of the complications arising from the nuclear-spin 
multiplicity, 'the quantum weight of rotational states of a 
molecule is worth special discussion. In general q r for any 
state = nuclear spin degeneracy x rotational degeneracy. For 
a rotational quantum number J , the rotational degeneracy is 
always (2J +1). The nuclear-spin degeneracy for molecules 
composed of dissimilar atoms is always 

^nucl " (2& s -f~l)(24'+l), 

where i s , i s > are the nuclear spins of the atoms in units hjZrr. 

As has already been pointed out, the nuclear-spin degen- 
eracy contributes a term to the T.N.S. functions derived 
from Z , which is not present in the H.T. functions. It is 
readily seen that in a reaction 

AB+CD ^ AC+BD, 

Z nucl AB = (2» A +l)(2i B +l), ^nuciOD = (2<o+l)(2i D +l), 
^nuciAC = (2i A +l)(2«o+l), Z nmi sl BD = (2iB+l)(2i D +l), 
^ciAB Z nucl CD 
Z nucl BB Z nncl AC ' 

and has no effect on A F = ET 2 hi G, so that for dis- 
similar molecules the nuclear-spin degeneracy term cancels 
in any chemical process. 

The NO molecule. 

For diato mi c molecules whose ground state is the rota- 
tional partition function is quite straightforward, viz. 

z mt = 
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When the molecule has more than one electronic state com- 
parable with JcT, two further points have to be noted in con- 
nexion with rotational quantum weights. Thus NO has two 
electronic states, 2 n^, 2 n&, separated by only Av . — ' 120 cm.™ 1 
The first point is that the smallest value of J = J for 2 n^ 
and J = f for 2 U%. The second is that, except for molecules 
in S states, each rotational level is subject to A-type doubling 
into two slightly separated levels whose energy is sufficiently 
close merely to double the quantum weight (2J-f-l). Since 
for oxygen i s = 0 and for nitrogen i s = 1, the nuclear-spin 
partition function is 1.3 — 3 — and the complete rota- 

tional partition function allowing for A-type doubling is 

^rot=3.2 2 (2J+1)« j' kT + 3.2 2 (2J+l)e-q**. 

J-*,#,!... b 

Rotational function for diatomic molecules made up of similar 
atoms . 

When the molecule is composed of two similar atoms, the 
thermodynamic complication may arise, as for ortho- and 
para-hydrogen at low temperatures, that transitions between 
odd and even rotational states do not occur. This possibility 
has already been discussed (cf. p. 102). 

When the ground state of the molecule is the following 
points have to be regarded when constructing the rotational 
partition function. 

If the nuclei have an even number of units of nuclear spin 
( \hj 27 r) they obey Einstein-Bose statistics, and the nuclear-spin 
multiplicity of rotational states with J even is (^+l)(2^-j-l), 
and with J odd is i s (2^+l)- 

An important example is deuterium, which has 

i $ = 2 units — 1, 

and the rotational partition function 

£rotn,= 6 2 (2J+l)e-^+«+3 T (2J+l)e~°W+V 

J even. Jodd 
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Hydrogen, on the other hand, has i s = 1 unit, and the 
nuclear-spin multiplicity of terms with J even is i s (2i s + 1) 
and with J odd (t s +l)(2^+l), i.e. the reverse of the above. 
The rotational partition function for hydrogen is thus 

£rotH 3 = 3 2 (2J+l)e~*W*»+ 2 (2/+ l)e-W+D. 

J odd J oven 

The molecule 0 16 0 16 . 

Since for oxygen i s — 0, the nuclei obey the Einstein-Bose 
statistics, and the alternate rotational states with J even 
have nuclear multiplicity 0, i.e. are missing. The ground 
state of the oxygen molecule is a 3 Sy state; the separation 
of the three energy levels of the triplet is negligible at ordinary 
temperatures, but the only rotational levels permitted for 
the three levels of the triplet are 

J = k+1, F % : J = k, F 3 : J — k — 1, 

where only odd values of k are permitted, so that 
^rot = 2 (2/C+3 )e-**J**+ 2 (2 K +l)er*«J**+ 

k — 1 , 3 , 5 ... 

+ 2 ( 2/c — i ) e~ € F3 * !kT . 

The electronic levels of the 0 16 0 16 molecule proximate to the 
ground level only begin to affect the complete partition 
function in the neighbourhood of 1,000° K. and 2,500° KL. 

Eotational partition functions at temperatures where a <^T. 

When a <^T 3 the various summations in these partition 
functions for symmetrical molecules can be replaced by 
integrations, and since in general 

(2J-\-\)e~ aJ( - J+v> = J (2/+1) dJ e-o-JV+U 

0 

co 

= J e~ ax dx = 1/a, 


where 


x = J(J+ 1 ), 
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it will be seen that for symmetrical diatomic molecules the 
complete rotational partition function including nuclear-spin 
multiplicity is (2* s +l) 2 1 

"rot 


It must be noted that in a reaction where unsymmetrical 
diatomic molecules are produced from symmetrical mole- 
cules, such as 

H 2 +D 2 ^ 2HD (cf. p. 110), A 2 +B 2 ^ 2AB, 

a consideration of the T.N.S. available energies, 




— JX1 m - |^— ^vib "elect "trails J 

-RTln {2iA f 1)2 G' At , 


A Bt = -i?yin (2 -^± 1 l 2 ^„ 


but 2A ab = — RT ln(2i A -j- 1 ) 2 (2i B -j- 1 ) 2 G AB 
and A A = -RT ln[2x2] Y%,~ , 

shows that the nuclear-spin multiplicity always cancels, 
except for the symmetry factor 72 In 2 for symmetrical 
molecules, so that H.T. functions are obtained by omitting 
the nuclear-spin factors, but retaining the symmetry factor, 
as was done in the expression for the theoretical chemical 
constant of diatomic molecules (cf. p. 37). 

The only important cases where a is not small compared 
with T occur for hydrogen below 273° K. and for deuterium 
below 200° K., and have already been discussed. 


The OH radical. 

The ground level for the OH radical is an inverted 2 II term, 
i.e. energy of 2 II$ < 2 n*. The separation of each of these 
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into A and B levels is large, so that instead of multiplying 
each rotational partition function by 2, as for the A-type 
doubling of NO, it is better to take four different partition 
functions. The nuclear-spin multiplicity is due to i a = J 
for H and i B = 0 for O, i.e. is 2. 

Z I0t = 2 [ 2 (2 J+l)e- € J lkT 1 for all values of ej. 

Examples for other diatomic molecules could be quoted, 
but show no special features of interest. Eor polyatomic 
molecules the calculation becomes very difficult, since the 
spectra have by no means been completely resolved. An 
approximate solution is possible by treating the energy states 
as independent and the molecules as rigid, so as to separate 
the partition function into factors. The vibrational partition 
function is not always easy to evaluate, particularly for long 
molecules. Eor smaller molecules it is relatively unimportant, 
e.g. for N a O at 298° K. Z Tot — 496, Z^t, = I’l- 

Linear molecules. HCN, N 2 0, C0 2 , C 2 H 2 . These are treated 
as diatomic, and Z rot = 1/scr, where s is the symmetry 
factor, and integration is permissible since T is large and I 
is large. 

Spherical rotators. CH 4 , CC1 4 have I A = I B = I c , and 
when integration is permitted Z Iot = V7 t/sct*. It may be noted 
that s = 12 for these two molecules, (s is defined as the 
number of distinguishable permutations produced by rota- 
tion of the molecule or its parts.) 

Symmetrical top molecules . NH 3 , CHC1 3 , C 2 H 6 (double top) 
have I A = I B j Ci and when integration is permitted an 
approximate expression for Z TOt is 

^ _ Vtt 

For the double top Z TOt — 

Sa A a O 
B 


4360 



122 THERMODYNAMICS OF HYDROGEN 


Asymmetrical top molecules. I c . A convenient 

approximation for this case is to write 

I = V(44U z lot = - It — - — 

SN\cr A a s o c J 

Typical values of s are H 2 0 = 2, C 6 H 6 = 12. 

For more complex molecules the expression 

SIT 

has been suggested, but direct experimental verification is 
difficult. 

References to substances whose thermodynamic functions 
have been evaluated spectroscopically are listed in Annual 
Reports Chem. Soc. 32, 84 (1936). Amongst the equilibria 
not easily accessible by other measurements the following 
may be quoted: 

H 2 0 = £H 2 +0H. N 2 = 2N. NO = N+O. C+2H, = CBL 

*0+2H»»0 t H t . 2C+H 2 = C 2 H 2 . Br 2 = 2Br. I 2 = 21. 

P 2 = 2P. S 2 = 2S. SO = JS 2 +J0 2 . 

It may be noted that thermodynamic functions calculated 
spectroscopically agree very well with those evaluated in 
other ways, with one or two exceptions, such as NO and CO, 
already mentioned in connexion with the verification of the 
heat theorem, and with the exception of substances such as 
ice and ethane, where the discrepancy may be due to incom- 
pletely investigated spectroscopic energy levels. 
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THE THERMODYNAMIC FUNCTIONS OF BLACK- 
BODY SPACE 

If a cavity is made in a solid body in temperature equi- 
librium, a small amount of energy will always be present 
in the cavity even when it is completely freed from matter. 
This is due to the presence of electromagnetic radiations in 
the ether. When equilibrium has been reached, the distribu- 
tion amongst the possible wave-lengths in the ether merely 
depends on the temperature, and not on the nature of the 
surrounding bodies. Owing to the fact that a perfectly black 
body (with emissivity unity for all wave-lengths) would 
emit radiation with exactly the same energy distribution, 
the radiation in a cavity in temperature equilibrium with its 
surroundings is usually called black-body radiation. 

Since the energy content of the ether can vary like that of 
any material substance, its thermodynamic functions can 
be calculated, not in terms of unit mass, but unit volume. 
It may be noted that it is desirable to extend the analogy 
between matter and radiation as far as possible in view of the 
fact that at very short wave-lengths there is a possibility 
of the process 2 hv = e + +e_. 

If E is the energy content in volume V , electromagnetic 
theory shows that the equation of state of the ether is 

PV = IE. 

The relation between the total energy (£ per unit volume of 
a black-body enclosure and the temperature may be estab- 
lished by ordinary thermodynamical reasoning. 

The entropy change in a reversible isothermal expansion 
at constant pressure is 

ds= mv+m = l^+py r+ y is . 

R 2 
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THE THERMODYNAMIC FUNCTIONS 


Since dS is a perfect differential (cf. p. 3), the R.H.S. must 
also be a perfect differential, i.e. 

_a/(£+P\ _ a V 1 
0<E\ T ] 8V T T' 

Substituting P = (£/ 3, 


d 

(m 

l 

d&\3T) 

~ T’ 

4 (S dT 

4- 4 

1 

3 T 2 d<£ 

+ 3 T 

= 

±dT 

d<£ 


T 

~~¥’ 



and (£ = aT* (Stefan’s law), 

where a is a constant. 

The entropy is thus 

dS = iaT*dV+±aVT*dT 
= iad(T*V), 

8 = iaVT*+b, 
where b is a constant. 

Distribution of energies among the different wave-lengths in 
black-body radiation. 

This distribution can only be calculated from statistical 
theory. The problem is how to distribute the total energy (£ 
amongst the different wave-lengths so as to obtain the most 
probable distribution. 

The number of stationary vibrations in a volume V between 
wave-lengths A, A+dA is 

dZ A =z 87tA -4 dXx F,t 

or since vA = c, 

dZ v = ^^-dv per unit volume. 

t This number is twice that on p. 64 to allow for the possible planes 
of polarization of the electromagnetic waves. 
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The distribution may be calculated from the Einstein-Bose 
statistics, or more simply from Planck’s original assumption 
that each of these stationary oscillations of the ether can only 
be associated with energy 

= h, 2hv, nh, 

where n is an integer and h is Planck’s constant. At tempera- 
ture T the probability of a vibration of frequency v with 
energy hv, 2 hv, and so on, is given by the Boltzmann formula 

m — Rp-nhvIkT 
rnv 

The mean energy of a vibration of frequency v at tempera- 
ture T will be obtained from the partition function (cf. p. 61), 
Z = J e~ nhvlkT , 

hv 

= f 


and since there are dZ v such vibrations the quantity of 
energy of the black-body radiation contained between v and 
v+dv is 


' I. 


e v dZ v 


8ttAv 3 dv 



XII 

ADSORPTION EQUILIBRIA 

The adsorption equilibria of gases on solids are worth dis 
cussing both because they represent a very characteristic 
physico-chemical equilibrium, whose thermodynamical im- 
plications are therefore worth attention, and because they 
illustrate how equilibrium formulae may be obtained from 
the notion of dynamic equilibrium, and the velocities of 
sorption and desorption. 

The Langmuir adsorption isotherm . 

The only case of adsorption to be discussed here is that 
of the ‘monomolecular layer 5 of gas held at certain adsorbing 
points on the solid, which are assumed to be sufficiently far 
apart for the adsorbed molecules to have no influence on one 
another, leading to an adsorption potential </> independent 
of the number of molecules already absorbed. 

Dynamic adsorption equilibrium. Equilibrium will be estab- 
lished when the number of molecules evaporating in unit 
time equals the number of molecules condensing. If there are 
N spaces in all available for adsorption, and the fraction a 
is occupied, the number evaporating in unit time will be 

§ ivenb y aocNe-^T, 

where e~ El ^ T is included to allow for the fact that the mole- 
cules must possess energy E x above the average before they 
can escape, and a is a constant. 

The number of molecules condensing will be 
bp ( 1 — ol)N e- E *l kT , 

since the number approaching the surface in unit time is bp 
(b = ^](I/27rmJcT), where m is the mass), and the fraction of 
the surface uncovered is ( 1 — oc)N. The term e~~ E ^ kT is included 
to allow for the fact that the molecules may have to possess 
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some activation energy before they will adhere on collision 
with an empty space. 

Equating the number of molecules condensing with the 
lumber evaporating, 

Naoc = bp(l—oc)Ne-^~ E ^ kT = bp{l— <x)Ne-& kT . 
Solving for a, 

bpe~& kT 

no. of molecules adsorbed = Na = N — — r ttt^. 

a+bpe~w kT 

This formula has the characteristic property that when p is 
small, the number of molecules adsorbed is proportional to p 

= N^e-V^p, 

and when p is large the number adsorbed is 
N = constant, 

i.e. the solid shows the phenomenon of saturation. In actual 
fact the assumption about non-interference of the adsorbed 
molecules usually becomes quite unsound long before satura- 
tion is reached. 


Statistical calculation . 

The advantage of a statistical over a dynamical calcula- 
tion is that the magnitudes of all the constants occurring in 
the formula are at once determined. 

Eor the distribution of molecules between sorbate and gas, 
the calculation may be made as follows : 

The probability of finding any molecule in the solid 
sorbate is = 

where is the maximum number of molecules which 

could be adsorbed, and N a is the number actually adsorbed, 
so that N m ^ x — N a = the number of available spaces. The 
probability of finding it in the gas phase is 



ir (27rmkT)% 

gas ^3 > 
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N a ^ N m&x -N a ¥ 


SO that » o-<j>jkT 

N g V g (2imkT)i ' 

On rewriting this equation it is seen to give the same 
as before, namely, since V g jN g = kT’jp, 

N^e-tWbp 


where 


N„ 


b = 


l J rbpe~ ( l > l kT 

¥ 1 


{2-nmkTf kT‘ 


The comparison just made between the statistical 
dynamical methods of calculating equilibrium illustrates 
general features of the dynamical method of calculai 
any equilibrium. When the phenomena occurring are fa 
simple, the correct dependence of the equilibrium on pressi 
etc., is obtainable without difficulty and with much sim] 
theoretical foundations than in the statistical method.; 
is usually necessary to leave the velocity constants undel 
mined, however. 

The statistical method, on the other hand, not only dei 
mines the dependence of equilibrium on pressure, etc., 1 
also gives the ratio of the velocity constants of the proces 
in dynamic equilibrium. It is from this property that exb 
sions of thermodynamic theory to the velocity of reactic 
are being developed. 
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